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Statistical Investigation of the Fatigue Life of 
Deep-Groove Ball Bearings 

J. Lieblein and M. Zelen 

Fatigue is an important factor in determining the service life of ball bearings. Bearing 
manufacturers are therefore constantly engaged in fatigue-testing operations in order to 
obtain information relating fatigue life to load and other factors. Several of the larger manu- 
facturers have recently pooled their test data in a cooperative effort to set up uniform and 
standardized ball-bearing application formulas, which would benefit the many users of anti- 
friction bearings. These data were compiled by the American Standards Association, which 
subsequently requested that the National Bureau of Standards perform the necessary 
analyses. This paper summarizes the principal results of the analyses undertaken by the 
Bureau, and describes the statistical procedures used in the investigation. 

1. Introduction 
1.1. Statement of Problem 

The experience of ball-bearing manufacturers over many years has led (o the 4 acceptance 
of an equation of the form [15, p. 15, eq (53)] 1 

L=(c/py, (i) 

relating fatigue life L to load P when other factors are kept constant. In the above equation, 
Cis termed the "basic (dynamic) capacity/' and is defined [15, p. 48] as the 4 constant bearing 
load (in pounds) that 90 percent of a group of similar bearings can endure for one million 
revolutions under tin 4 given running conditions. 

The quantity C in eq (1) depends upon the characteristics of the 4 bearing type, as indicated 
in [15, p. 32, eq (120)]. When the expression cited is substituted in eq (1), the fatigue-life 
formula for ball bearings takes the form 

L J f c Z^Dp(i cos a) a * -y> 

The s t ymbols are defined as follows: 
Z= number of balls. 
D a =bal\ diameter in inches. 
i= number of rows. 
a= contact angle. 
P=bearing load in pounds. 

L= number of million revolutions that a specified percentage of bearings will fail to sur- 
vive 1 on account of fatigue causes. If the percentage is 10, then L = L l() , and is 
termed the rating life; if the percentage is 50, then L=L 5Q , the median life. 
Pj (ti, a 2 , a d ,f c are taken as unknown parameters whose values have to be estimated from 

given data. 
Since i=\ and a = ()° for deep-groove ball bearings, with which this paper is exclusively 
concerned, the life equation that will henceforth be considered takes the form 



j f t .Z'iDi np 



L=\ J -^- a ~ ■ (2a) 



Figures in brackets Indicate the literature references al the end of this paper. 
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The main goal of this investigation was to determine the "best values" of the unknown 
parameters in the life equation from the experimental data. One of the major problems was 
to determine the value of the exponent p, as there was disagreement within the ball-bearing 
industry whether an appropriate value for p was 3, 4, or some other value. 

1.2. Description of Data 

The data available for analysis consisted of sets of records summarizing endurance tests 
for deep-groove ball bearings. These tests were carried out over a period of years by four 
major ball-bearing companies. In the interest of trade anonymity, these companies will 
henceforth be designated by A, B, C, and D. Each endurance test consisted of a number of 
bearings of the same type (the number varying from test to test), which were tested simul- 
taneously under the same load and running conditions. Table 1 summarizes the number of 
test groups of data for each company. The data from company B were sufficiently extensive 
to permit a further breakdown into three bearing types, here denoted by B-1, B-2, and B-3. 





Table 1. 


Summary of ball-bearing data 






Company 




Number of test 
groups 


Total number of 

bearings in test 

groups 


A 


50 
148 

37 
94 
17 
12 
3 


1, 259 
3, 289 


B 


Type B-1 

Type B-2 

Type B-3 

C 






291 

109 


D _ 


Total (all 


companies) 




213 


4, 948 



The worksheets, summarizing the tests, recorded the number of millions of revolutions 
reached by each bearing in the test group before fatigue failure. Information was also given 
for those tests terminated before all bearings in the test group failed. In addition to the test 
results, the worksheets included information on the characteristics of the bearing type (e. g., 
values for Z, D a , i, a) and load P, as well as other items of descriptive and identifying informa- 
tion. A specimen worksheet is reproduced in appendix A. 

All necessary quantities for evaluating the unknown parameters in the life equation (2) 
were given directly on the worksheets except the fatigue life L. 2 This quantity can be estimated 
from the observed fatigue lives of individual bearings within a test group. As already noted, 
two concepts of fatigue life are used for L, namely, the rating life L 10 , and the median life L 50 . 
Separate analyses have been carried out with regard to each throughout. 

Appendix A summarizes the data taken from the original worksheets that were used in 
the statistical analysis. Also given are the computed values for i 10 , Z 50 , and the "Weibull 
slope" e (which relates to the dispersion of fatigue lives). The methods for obtaining these 
quantities from the bearing data are given in detail in appendix B. 

1.3. Assumptions for the Statistical Analyses 

All conclusions reached in this report, and all statistical analyses employed, are based upon 
the following principal assumptions : 

(a) The life formula (2) is the proper functional form for describing fatigue life in ball 
bearings. 

(b) Differences in the measured life of bearings classed as identical, tested at the same load, 
reflect only the inherent variability of fatigue life, and are free from systematic errors that may 
arise from different test conditions, materials, manufacturing methods, etc. 

2 Certain estimates for L\o and Z50 had been entered on the worksheets for many of the tests. However, these were not regarded as part of the 
data submitted for analysis. 

274 



(c) AH the bearings in a test group can be regarded as a random sample from a homogeneous 
population of ball bearings. 

(d) The probability distribution of the number of revolutions to fatigue failure is of the 
same form for each test group, although its parameters may differ from group to group. 

(e) This fatigue-life distribution is of the type known as the "Weibull distribution/' 
Tli" purely statistical assumptions, (e) to (e), served as the basis for the determination 

of L 10 , L 50 , and e for each test group. Assumption (e), however, is not involved in the methods 
used to evaluate the parameters in the 4 life formula (2) from given values of i l0 or L m . A 
different assumed form for the distribution of fatigue life might give somewhat different values 
for Li and L m , but the same methods could then be used to evaluate the unknown parameters 
in the life formula (2). 

Other assumptions of a more technical nature were necessary in the course of the analyses. 
These are discussed in appendixes B and C. 

As in all cases where inferences are made from given data, the conclusions reached here 
pertain only to the population from which the given data can be regarded as constituting a 
random sample. 

2. Outline of Statistical Analyses 

The statistical analyses were divided into two phases. The first phase considered the 
problem of finding estimates of L in , L 50 , and the Weibull slope e from the given test data; the 1 
second phase used these estimates of L 10 and L m to evaluate the unknown values of the param- 
eters in the life formula. 

2.1. Estimation of L !0 and L m 

The quantity L depends upon the existence of an underlying probability distribution of 
bearing lives. Selection of a distribution or population is equivalent to specifying the proba- 
bility that a bearing selected at random from such a population will survive any given number 
of revolutions, />, or, conversely, that if c is a specified probability, then L is the life period 
that will be survived with this probability, e. g., 



Probability {life>Z}=c= 



.90 for L=L 10 . 
.50 for L = L 50 . 



Accordingly, any L, such as L U) or Z 50 , must be obtained by estimating a characteristic of 
the assumed distribution. For reasons described in appendix B, the distribution characterizing 
ball-bearing fatigue life was taken to be the Weibull distribution. In brief, this distribution 
can be derived by assuming a "weakest-link" concept of fatigue strength. In addition, the 
suitability of the Weibull distribution for fatigue life has been verified in many cases by empirical 
plotting of data. 

One method of estimating L l0 or L 50 makes use of special probability-plotting paper so 
designed that a theoretical Weibull distribution plots as a straight line, and treats the problem 
as one of straight-line fitting by conventional least squares procedures. However, t he procedure 
usually followed does not take into full account the number of bearings that remain intact when 
tests are incomplete, nor the interdependence of successive points. Because of these and other 
limitations, it was decided to use an alternative approach in the estimation of L 10 and L 50 for 
each test group (see appendix B). 

To this mid, a method was developed that takes into account explicitly the number of 
bearings remaining intact at the termination of a test, and that also possesses several other 
advantages. This method makes use of certain specially determined linear functions of the 
observed failure times (in logarithms), x t) arranged in order of size. These functions have the 
general form 

r=&*,. (3) 

i=i 
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As the method makes intimate use of the ordered arrangement of the data, it is termed an 
"order statistics" method. 

The coefficients Cj in eq (3) allow great flexibility. They have been determined in such 
a manner that the method will have certain desirable objective characteristics, e. g., freedom 
from systematic error and a minimum standard error. 

2.2. Evaluation of the Parameters in the Life Formula 

Once the estimates for i 10 and L 50 are obtained, it is possible to evaluate the exponent p 
in the life formula. However, in order to make the most efficient use of the given data, it is 
necessary also to estimate the other parameters, f c , a iy and a 2 . 

The methods for estimating the values of L l0 and Z 50 for each test group actually yield 
results for In L l0 and In L 50 . Thus, taking logarithms 3 of the life equation (2 a) gives 

In L= (p Jn/ C )+ (pa,) In Z+ (pa 2 ) In D a -p In P. (4) 

This equation can be written more simply as 

Y = bo+b 1 x l + b 2 x 2 + b 3 i 3 , (5) 

where 

F=m L (for either L i0 or Z 50 ), 



&o=pln f e , b 1 =pa 1 , b 2 =pa 2 , h=—p, 
^!=ln Z, x 2 =\nD a , ;r 3 =lnP. 



(6) 



The quantities x u x 2 , and x z depend on the characteristics of the bearing type and test 
conditions, and can be regarded as known exactly. On the other hand, the variable Y, which 
depends on the outcome of the bearing tests, is subject to considerable dispersion. Thus, 
estimates can be found for the parameters 6 , b u b 2j and 6 3 , using standard least squares methods 
based on minimizing the sums of squared deviations in the y direction. These methods are 
discussed in detail in appendix C 

After the parameters 6 , b u b 2 , and 6 3 are estimated, values for a , a u a 2j and p can be 
found from the relations 

i , b bi 

ao=ln; c =— t- ai=— } 



h h 



h 



a 2 =—^ p=-h 

It is clear that the values for a , a u and a 2 depend on the value of p. 

The estimates for p and the a's are subject to some uncertainties because they are based 
on test results, which themselves are subject to considerable variability. Hence with every 
value of p and of the a's calculated from the life data, there is given also an interval of un- 
certainty to indicate its precision. These intervals are " 95-percent confidence limits." 4 

A large interval of uncertainty associated with an estimate indicates poor precision; a 
small interval of uncertainty is evidence of high precision. These intervals of uncertainty 
not only reflect the inherent variability of the test data, but are also affected by (a) how well 
the life equation (2a) is the proper functional form for bearing life, and (b) the suitability of 
the data (including the number of test groups) for estimating the parameters in the life formula. 

Further technical details concerning the evaluation of the parameters in the life formula 
are given in appendix C. 

3 Natural logarithms to the base e are used throughout. 

4 Briefly, confidence intervals describe the compatibility of the observations with an unknown parameter estimated from them; 95-percent 
confidence limits are limits such that on the average, in repeated applications of the same procedure, 95 percent of intervals so calculated will 
contain the unknown true value of the parameter. The confidence limits associated with p are symmetric. However, the confidence limits asso- 
ciated with the a's are asymmetric because of the dependence of the a's on p. 
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3. Summary of Analyses 



3.1. Evaluation of Parameter p 

The statistical analysis based on all deep-groove ball-bearing data from companies A, B, 
and C 5 yielded the final values for p shown in tabic 2. The separate values for each of the 
three companies arc given in table 3. The intervals of uncertainty specified by ± quantities 
in those tables refer to intervals within which, with reasonable assurance, the true value of 
the parameter is located. The fact that all of the intervals of uncertainty exhibit considerable 
overlap shows that the data are consistent with the supposition that all three companies have 
a common value of p for deep-groove bearings. The fact that all the intervals include .'> 
indicates that all the estimates of p are consistent with the practice of taking p=3. More- 
over, the value of p for L w was not significantly different from that for L 50 . 

Table 2. Final over-all valves of p for deep-groove 
bearings 



V- 


LlQ 


L 50 


2. 87 ±0. 35 


2. 80 1 0. 31 



Table 3. Individual estimates of p for deep-groove bearings by company 



( Jompany 


Number 
of test 
groups 


L/IQ 


/>. r )0 


A 

B._ 
C__ 


50 

148 
12 


3. 00 ± 0. 64 

2. 75 ± . 48 

3. 12 ± .88 


3. 05 1 0. 60 
2. 62 d 0. 40 

2. 88 ± 1. 02 



The values given for p are based on analyses of all deep-groove ball-bearing data, irrespec- 
tive of bearing type. Hence, the parameter estimates represent "omnibus" values. In order 
to investigate the dependence of the exponent p on bearing type, (he data from company B, 
which was made up of three bearing types, were analyzed separately. The results for the 
exponent p are shown in table 4. These results are all compatible with the value j) = 'S. 

Table 4. Value of p by bearing type for company B 



Type 


Number 
of test 
groups 


Value of p 


L\o 


Lo0 


B-l...... 

B-2 

B-3____ 

Total 


37 

94 
17 


3. 36 ±0. 68 
2. 65 ±0. 91 
1. 89 ±1. 28 


3. 23 ±0. 47 
2. 13 4 0. 79 

2. 82 j 1. 10 


148 











5 The data furnished by company I) were too few to be Included in the analysis. 
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3.2. Evaluation of Parameters f c , a u and a 2 

The computations that give estimates for the exponent p also yield estimates for the quan- 
tities In f c , a u and a 2 . From the relations (6) it is clear that the values for these parameters 
depend on the value for p. Thus, associated with every value of p will be corresponding- 
values for In f c , a u and a 2 . Table 5 summarizes these parameter estimates associated with 
the final values of p. The estimates for a =ln/ c , rather than f C} are given here, because this is 
the parameter that arises naturally in the life formula (cf. eq 4). 

The analyses conducted separately for each company resulted in other values than those 
in the previous paragraph for a , a u a 2 . These results are summarized in table 6. They show 
excellent agreement with the results in table 5, even though the values for p are somewhat 
different. 







Table 5. Final vain 


es of a , «i, a 2 for L\ and L- n0 






Company 


V 


a 


Interval of 
uncertainty 


a x 


Interval of 
uncertainty 


a 2 


Interval of 
uncertainty 


£ao 


A 

B 

C 


2. 87 
2. 87 
2. 87 


9. 02 

8. 55 

9. 56 


(7.31, 10.79) 
(7.98, 9.14) 
(6.85, 12.42) 


0. 380 

. 670 

-. 174 


(-0.454, 1.201) 
( 0. 418, 0. 920) 
(-1.750, 1.352) 


1. 72 
1. 81 
1. 37 


(1.51, 1.92) 
(1.70, 1.92) 
(0. 09, 2. 67) 


^50 


A 

B 

C 


2. 80 
2. 80 
2. 80 


10. 36 
9. 05 
9. 05 


(8.81, 11.98) 
(8.54, 9.60) 
(6.61, 11.58) 


0. 015 
. 695 
. 475 


(-0.741, 0.751) 
( . 470, . 920) 
(- . 921, 1. 847) 


1. 69 
1. 91 
1. 76 


(1.50, 1.88) 
(1.81, 2.01) 
(0. 60, 2. 93) 



Table 6. Values of a , «i, «2 for L i0 and L 50 , based on independent analyses for each company 



Company 



a 



Interval of 
uncertainty 



Interval of 
uncertainty 



«2 



Interval of 
uncertaintv 



A 
B 
C 

A 
B 

C 



3. 00 


8. 97 


2. 75 


8. 59 


3. 12 


9. 21 



(7.18, 10.90) 
(7.99, 9.24) 
(7.29, 11.84) 






390 




666 


— 


041 



(-0.507, 1.249) 

( . 398, 0. 928) 
(-1.326, .992) 



1. 


73 


1. 


80 


1. 


36 



(1.50, 1.94) 

(1.67, 1.92) 
(0. 49, 2. 30) 



3. 05 


10. 13 


2. 62 


9. 15 


2. 88 


8. 93 



(8.48, 12.00) 
(8. 61, 9. 76) 
(6.58, 12.39) 



0. 072 
. 690 
. 510 



(-0.768, 0.855) 
( . 456, . 922) 
(-1.055, 1.810) 



1. 


71 


1. 


90 


1. 


75 



(1.50, 1.91) 
(1.79, 2.00) 
(0. 66, 3. 05) 



Similarly, the values for a , a u and a 2 , arising from separate analyses made on the three 
types of bearings from company B, resulted in still other estimates for these parameters. Table 
7 summarizes these estimates. These estimates are less precise than the corresponding omnibus 
values given for company B in table 6. This is a consequence of the fact that within a bearing 
type, the quantities Z and D a hardly vary at all. This condition makes the data unsuitable 
for estimating the associated unknown parameters, a , a u and a 2 . 
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Table 7. Values of a , «i, «2 .for Lio cmd L 50 6?y bearing type (company B) 



Type 


a 


Interval of 
uncertainty 


«! 


Interval of 
uncertainty 


a 2 


Interval of 
uncertainty 


Lio 


B-1 

B-2 

B-3____ 


7. 25 
7. 34 
2.50 


( 5.21, 9.39) 
( 5. 54, 9. 33) 
(-7.95, 16.04) 


1. 07 
1. 21 
3. 70 


( 0.35, 1.76) 
( . 43, 2. 00) 
(-2.25, 9.06) 


1. 68 
1. 69 
1. 27 


(1.27, 2.08) 
(1.38, 1.93) 
(0.30, 1.65) 


jL/50 


B-1 

B-2 

B-3 


7. 39 
9. 00 
1. 03 


( 5. 92, 8. 92) 
( 7. 10, 11. 67) 
(-4.34, 6.53) 


1. 23 
0. 87 
4. 50 


( 0. 73, 1. 73) 
(-.05, 1.68) 

( 1. 97, 7. 19) 


1. 79 
1. 77 
1. 48 


(1.50, 2.08) 
(1.46, 2.03) 
(1.21, 1.70) 



3.3. Redetermination of the Estimates for f c 

The uncertainty intervals associated with estimates for the parameter ao=ln/ c are quite 
large. This is primarily because the uncertainty associated with the estimate of a also de- 
pends on how well the other parameters, a u a 2 , and />, are estimated. Another way to evaluate 
a , which may result in smaller intervals of uncertainty, is to assume a priori values for a u a 2l 
and p, and then determine the estimate for a {) . This procedure was followed by using the 
widely accepted values for the parameters given in [15], namely, ai=2/3, a 2 =1.8, p=3. 

However, if on such a calculation the values assumed for the parameters a u a 2 , and p are 
not compatible with the given data, then values of a () (or f c ) so calculated will not be correct 
determinations for these data. Accordingly, an analysis was made to determine whether the 
parameter values in [15] were compatible with the given data. 

This analysis showed that these parameter values are compatible with the data, with re- 
spect to all individual companies for rating life L 10 , but not for median life /> 50 . (Company 
A was the only company for which the parameter values are suitable for median life.) A fur- 
ther analysis, by bearing type tor company B, showed that the above parameter values are 
not suitable for the rating life L x[) with respect to B-3-type bearings. 

In the light of this last analysis, redetermined values of &o, taking ai = 2/3, a 2 =1.8, and 
p=3, are only strictly valid with respect to company A, company B (B-1, B-2), and company 
C for rating life X ]0 . These values are summarized in table 8. For convenience, these new 
estimates are given for/ c =ln _1 ao=exp a . 

Table 8. Values for f c assuming ai = 2/3, a 2 = 1.8, p = 3.0 for L 10 



Company 


Number of 
test groups 


fc 


Interval of uncer- 
tainty 


A_ 


50 
148 
37 
94 
17 
12 
3 


4,538 
4, 925 
4, 709 
5,033 

3," 294 

4, 639 


(4,273, 4,817) 
(4,750, 5,105) 
(4,403, 5,034) 
(4, 885, 5, 187) 


B (over-all value) 

B-] 


B-2__ 


B-3*_ 


C_ 


(3,029, 3,583) 

(3,478, 6,187) 


D 





^Assumed values of parameters a\, a 2 , and p not compatible with test results for bearings of this series. 
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4. Appendix A. Summary of Original Data 

This appendix summarizes in tabular form the worksheets submitted by the American 
Standards Association Subcommittee to the National Bureau of Standards for statistical 
analysis. Separate tables are presented for deep-groove data from companies A, B, C, and D. 
These four tables (A-l to A-4) are followed by table A-5, which gives a synopsis of the number 
of test groups and the number of bearings for each company. 

Tables A-l to A-4 give the size of test group, the values for quantities P, Z, D a , and the 
estimates 6 for i 10 , i 50 , and the "Weibull slope" e. All of these variables are directly observed 
or specified quantities except for the estimates i 10 , i 50 , and e. These last three quantities are 
based on statistical calculations that made use of the results of individual endurance tests. 
These calculations are explained in appendix B. 

The original data, as submitted, contained a few cases where companies tested bearings 
manufactured by other companies. Such test groups are not included in the summary tables, 
as these results confound differences in testing with differences in manufacturing. Therefore 
these test results were not used in any of the analyses. Thus, table A-3, for company C, 
omits 4 tests performed on other manufacturers' bearings; table A-4, for company D, omits 3 
tests. 

The five tables described above are followed by a specimen worksheet " with identifying 
information removed. A sample of Weibull-f unction coordinate paper is also included. This 
coordinate paper had been used for graphing the results of all the individual endurance tests 
and these graphs had accompanied the worksheets submitted to the Statistical Engineering 
Laboratory. 



6 The estimates for Zio and £50 are given in millions of revolutions for all companies except company D. The life estimates for D are shown in 
hours, the same units in which the original endurance data were given. 

7 Bearings marked "Omitted" were completely eliminated from consideration, as company representatives explained that these were non- 
fatigue failures and should not be regarded as part of the test group. As a result, the test group in the case of the specimen sheet shown was taken 
to consist of 23 bearings rather than the original number of 25. This type of situation appeared rather infrequently, however. 
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Table A-l. Summary ball-bearing data for company A, with computed values for L i0j L 50 , and Weibidl slope e 



Record 


Year 


Number 




Z 


D a 






e 


No. 


of 


in test 


Load 


Number 


Ball diam. 


L 


L 5o 


Weibull 




test 


group 




of balls 




10 


slope 








lb 




in. 








1- 1 


1936 


2k 


1*22*0 


8 


11/16 


19.2 


Sk.5 


1.27 


1- 2 


1937 


20 


U2k0 


8 


11/16 


26.2 


7l*.2 


1.81 


1- 3 


1937 


111 


U2U0 


8 


11/16 


11.1 


68.1 


1.0)4 


1- k 


1937 


19 


U2U0 


8 


11/16 


11.8 


66.8 


1.09 


1- 5 


1937 


18 


U2U0 


8 


11/16 


13.5 


79.1* 


1.06 


1- 6 


1938 


21 


2530 


9 


1/2 


5.80 


25.7 


1.27 


1- 7 


1938 


28 


1*21*0 


8 


11/16 


18.3 


kk.l 


2.10 


1- 8 


1938 


27 


U2U0 


8 


11/16 


5.62 


73.2 


.73 


1- 9 


19U0 


20 


U2U0 


8 


11/16 


15.8 


82.7 


1.12* 


1-10 


19k0 


22 


U2U0 


8 


11/16 


8.70 


m.6 


1.20 


1-11 


19h0 


19 


U2U0 


8 


11/16 


11.6 


160 


.72 


1-12 


19k0 


15 


19i*0 


9 


1/16 


20.6 


71.U 


1.52 


1-13 


19h0 


15 


19k0 


9 


7/16 


lli.5 


88.2 


1.0I* 


l-ll* 


19k0 


15 


2536 


9 


1/2 


12.1 


33.1 


1.87 


1-15 


19h0 


li* 


2536 


9 


1/2 


15.1 


1*6.1* 


1.67 


1-16 


19h0 


15 


2536 


9 


1/2 


ll*.0 


1*3.6 


1.66 


1-17 


19)40 


1U 


2536 


9 


1/2 


19.3 


51.8 


1.91 


1-18 


19h0 


26 


U2U0 


8 


11A6 


U6.2 


no 


2.17 


1-19 


19h0 


U* 


U2U0 


8 


11A6 


30.0 


88.2 


1.71* 


1-20 


191*2 


20 


U2U0 


8 


11/16 


21.1 


57.i1 


1.89 


1-21 


19h2 


20 


U2U0 


8 


11 A6 


17.3 


1*5.7 


1.91* 


1-22 


19h2 


37 


U2U0 


8 


11/16 


37.5 


118 


1.61* 


1-23 


19h2 


36 


1*21*0 


8 


11A6 


20.3 


77.1 


I.I4I 


L.2U 


19h2 


32 


1*21*0 


8 


11/16 


!*.03 


1*2.5 


.80 


1-25 


19l4l 


28 


2$kk 


8 


17/32 


8.38 


81*.7 


.81 


1-26 


191*3 


23 


3975 


8 


19/32 


1.79 


13.5 


.93 


1-27 


19h2 


30 


UUoo 


10 


5/8 


11.7 


1*5.1 


1.39 


1-28 


19k2 


31 


6920 


8 


7/8 


2*.i5 


15.8 


1.1*1 


1-29 


19h3 


30 


990 


9 


^/16 


7.23 


1*1.0 


1.09 


1-30 


19h3 


30 


1509 


7 


1/16 


22.9 


110 


1.20 


1-31 


19k3 


30 


932 


7 


11/32 


9.51* 


31.6 


1.57 


1-32 


192*2* 


26 


3180 


8 


19/32 


6.28 


23.0 


1.1*5 


1-33 


19hk 


29 


3180 


8 


19/32 


l*.8l 


21.2 


1.27 


1-31* 


19hh 


33 


86I1O 


10 


7/8 


1*.17 


12.8 


1.68 


1-35 


192*1* 


26 


11*080 


8 


1-1/1* 


5.1*2 


31.6 


1.07 


1-36 


1951 


28 


191*0 


9 


7A6 


7.1*7 


1*9.5 


1.00 


1-37 


1951 


3k 


2330 


9 


7A6 


lu 80 


21.3 


1.26 


1-38 


1951 


27 


1550 


9 


7A6 


ll*.8 


78.1* 


1.13 


1-39 


1951 


29 


1165 


9 


7A6 


8U.9 


1*60 


l.ll 


1-2*0 


1951 


27 


2910 


9 


7A6 


3.1*0 


16.5 


1.19 


i-ia 


1951 


27 


3880 


9 


7A6 


1.21* 


3.23 


1.97 


14*2 


1951 


26 


776 


9 


7A6 


22*1 


951 


1.37 


1-2*3 


1951 


30 


19750 


8 


1-3A 


3.01 


12.6 


1.31 


1-1*2* 




30 


2112 


8 


11A6 


89.1 


1*86 


1.11 


1-1*5 


a 


30 


1*221* 


8 


11/16 


15.2 


102* 


.98 


1-2*6 


0) 

1 


30 


8UU8 


8 


11A6 


2.0U 


10.2 


1.17 


1-1*7 


W) 


30 


2112 


8 


5/8 


51.0 


376 


.9k 


1-1*8 


-p 

3 


30 


1*221* 


8 


5/8 


5.26 


58.8 


.78 


1-1*9 


*£* 


30 


81*1*8 


8 


5/8 


.883 


k.9k 


1.09 


1-50 


191*2* 


30 


U22U 


8 


11A6 


11*. 8 


57.U 


1.39 
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Table A-2. Summary ball-bearing data for company B, with computed values for Lm, L 50 , and Weibull slope e 



Record 


Year 


Number 




z 


D a 






e 


No. 


of 
test 


in test 

group 


Load 


Number 
of balls 


Ball diam. 


\o 


L 5o 


Weibull 
slope 








lb 




in. 








2- 1 


191*0 


19 


570 


10 


3A6 


6.68 


13.1* 


2.72 


2- 2 


191*1* 


20 


570 


9 


5/16 


29.8 


70.0 


2.22 


2- 3 


191*6 


23 


580 


9 


1/1* 


16.3 


55.1 


1.55 


2- 1* 


191*6 


23 


580 


9 


1A 


28.5 


69.2 


2.13 


2- 5 


191*7 


23 


580 


9 


iA 


16.1* 


1*9.3 


1.71 


2- 6 


191*3 


10 


665 


9 


iA 


10.3 


1*0.1 


1.50 


2- 7 


191*1* 


10 


665 


9 


1/1* 


25.7 


1*6.1* 


3.19 


2- 8 


191*2 


19 


580 


10 


1/1* 


9.55 


39.6 


1.32 


2- 9 


191*6 


33 


620 


10 


iA 


11.9 


62.1 


1.51 


2-10 


191*7 


15 


620 


10 


iA 


19.9 


73.2 


1.1*5 


2-11 


191*7 


31 


620 


10 


1/1* 


12.9 


50.1* 


1.39 


2-12 


191*1* 


19 


625 


10 


1/1* 


19.3 


1*6.2 


2.18 


2-13 


191*1 


17 


720 


10 


1/1* 


11.1 


23.3 


2.51* 


2-1U 


191*6 


60 


980 


11 


9/32 


15.7 


1*3.5 


1.85 


2-15 


191*7 


32 


980 


11 


9/32 


11.2 


38.1 


1.5U 


2-16 


1950 


k9 


600 


11 


5A6 


1*17 


809 


2.85 


2-17 


191*9 


60 


600 


11 


5/16 


216 


709 


1.58 


2-18 


191*3 


20 


900 


11 


5/16 


35.6 


100 


1.82 


2-19 


191*6 


67 


1220 


11 


5/16 


12.0 


1*2.2 


1.50 


2-20 


191*7 


31* 


1220 


11 


5/16 


8.53 


1*6.6 


1.11 


2-21 


191*0 


20 


1370 


11 


5/16 


6.77 


18.9 


1.85 


2-22 


1950 


60 


11*15 


11 


5/16 


13.5 


h6.5 


1.53 


2-23 


1950 


60 


221*3 


11 


5/16 


2.32 


8.06 


1.51 


2-2l* 


191*2 


20 


720 


12 


5/16 


36.7 


11*1 


1.1*0 


2-25 


191*6 


55 


1300 


12 


5/16 


19.0 


57.2 


1.71 


2-26 


191*7 


30 


1300 


12 


5/16 


19.5 


60.6 


1.67 


2-27 


19hh 


20 


1650 


11* 


11/32 


17.0 


7l*.l* 


1.37 


2-28 


191*6 


$9 


1760 


11* 


11/32 


20.9 


53.7 


2.00 


2-29 


191*7 


31* 


1760 


11* 


11/32 


9.56 


1*0.7 


1.30 


2-30 


191*0 


20 


2010 


13 


13/32 


5.1*9 


33.3 


1.05 


2-31 


191*0 


9 


2010 


13 


13/32 


1.39 


l*l*.o 


.51* 


2-32 


191*1* 


19 


211*0 


13 


13/32 


9.80 


82.7 


.88 


2-33 


191*3 


11 


2630 


15 


13/32 


5.19 


51*.9 


.80 


2-31* 


191*2 


12 


5900 


11* 


19/32 


6.36 


17.5 


1.86 


2-35 


191*7 


19 


5900 


11* 


19/32 


3.68 


22.1 


1.05 


2-36 


191*7 


20 


8070 


11* 


23/32 


8.31* 


23.6 


1.81 


2-37 


191*2 


12 


8075 


11* 


23/32 


6.78 


36.1* 


1.12 


2-38 


1938 


10 


565 


9 


.210 


9.27 


18.1* 


2.75 


2-39 


191*0 


23 


720 


8 


9/32 


18.2 


56.9 


1.66 


2-UO 


191*0 


2U 


720 


8 


9/32 


22.8 


56.2 


2.09 


2-1*1 


191*0 


25 


720 


8 


9/32 


3.99 


15.6 


1.38 


2-1*2 


191*0 


21 


720 


8 


9/32 


9.07 


29.U 


1.60 


2-1*3 


191*0 


25 


720 


8 


9/32 


7.11i 


28.5 


1.36 


2-1*1* 


191*0 


25 


720 


8 


9/32 


12.5 


26.1* 


2.51 


2-1*5 


191*1 


25 


720 


8 


9/32 


18.8 


1*8.7 


1.98 


2-1*6 


19U1 


23 


720 


8 


9/32 


21.5 


53.2 


2,08 


2-1*7 


191*7 


33 


860 


8 


5/16 


17.1 


59.0 


1.52 


2-1*8 


191*8 


8 


860 


8 


5/16 


15.2 


87.6 


1.08 


2-U9 


191*3 


20 


900 


9 


5/16 


30.1 


92.3 


1.68 


2-50 


191*1* 


18 


900 


9 


5/16 


15.0 


1*7.6 


1.63 
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Table A-2. Summary ball-bearing data for company B, with computed values for L l0l L 50 , and Weibull slope e — Continued 



Record 


Year 


Number 




Z 


\ 






e 


No. 


of 


in test 


Load 


Number 


Ball diam. 


*lo 


L 5o 


Weibull 




test 


group 




of balls 






slope 








lb 




in. 








2-51 


191*5 


27 


91*0 


9 


5A6 


17.5 


52.8 


1.71 


2-52 


19l*7 


3h 


9h0 


9 


5/16 


li*.l* 


65.6 


1.21* 


2-53 


1938 


10 


II80 


9 


5A6 


8.76 


22.1 


2.01* 


2-51* 


191*5 


30 


1580 


9 


3/8 


12.1 


1*3.3 


1.1*7 


2-55 


191*7 


33 


1580 


9 


3/8 


17.2 


6h.6 


1.1*2 


2-56 


191*8 


8 


1580 


9 


3/8 


10.7 


3U.6 


1.61 


2-57 


19^ 


31 


2160 


9 


7A6 


10.9 


37.6 


1.52 


2-58 


191*7 


30 


2160 


9 


7A6 


12.7 


53.7 


1.30 


2-59 


1938 


9 


2200 


9 


7A6 


3.73 


1*3.5 


.77 


2-60 


191*7 


30 


21*80 


9 


7A6 


16.6 


78.3 


1.21 


2-61 


1950 


1*0 


1340 


9 


15/32 


180 


275 


1*.1*1* 


2-62 


1937 


19 


1660 


10 


7A6 


85.2 


23U 


1.86 


2-63 


19l|l 


19 


1700 


9 


15/32 


57.1 


230 


1.35 


2-61* 


1939 


21* 


21*80 


9 


15/32 


15.7 


55.8 


1.1*8 


2-65 


1939 


25 


21*80 


9 


15/32 


27.1 


97.8 


1.1*7 


2-66 


1939 


23 


21*80 


9 


15/32 


21.7 


122 


1.09 


2-67 


1939 


28 


21*80 


9 


15/32 


13.2 


1*2.3 


1.62 


2-68 


1939 


28 


21*80 


9 


15/32 


35.8 


11*5 


1.35 


2-69 


P-939 


20 


21*80 


9 


15/32 


12.7 


3U.7 


1.87 


2-70 


19hh 


20 


21*80 


9 


15/32 


10.1 


27.8 


1.87 


2-71 


19l*5 


20 


21*80 


9 


15/32 


8.83 


3U.3 


1.39 


2-72 


1938 


10 


21*80 


9 


15/32 


16.5 


60.3 


1.1*5 


2-73 


19142 


11 


21*80 


9 


15/32 


17.9 


65.8 


1.1*5 


2-7li 


191*3 


10 


21*80 


9 


15/32 


15.7 


63.1 


1.35 


2-75 


191*3 


20 


2U80 


9 


15/32 


10.8 


1*2.1 


1.38 


2-76 


19U* 


18 


21*80 


9 


15/32 


ll*.2 


39.9 


1.83 


2-77 


191*1* 


18 


2) 4 80 


9 


15/32 


19.0 


67.8 


1.1*8 


2-78 


191*1* 


18 


21*80 


9 


15732 


16.3 


57.7 


1.1*9 


2-79 


19W* 


20 


21*80 


9 


15/32 


2.93 


18.0 


1.01* 


2-80 


191*1* 


20 


21*80 


9 


15/32 


5.69 


25.1* 


1.26 


2-81 


19hh 


28 


21*80 


9 


15/32 


9.51* 


39.9 


1.32 


2-82 


±9hh 


22 


21*80 


9 


15/32 


12.6 


55.7 


1.27 


2-83 


I9hh 


23 


21*80 


9 


15/32 


5.10 


37.5 


.9h 


2-81* 


191*1* 


18 


21*80 


9 


15/32 


16.0 


53.7 


1.56 


2-85 


191*U 


20 


21*80 


9 


15/32 


1.98 


22.1 


.78 


2-86 


191*5 


20 


21*80 


9 


15/32 


5.65 


28.8 


1.16 


2-87 


191*5 


20 


21*80 


9 


15/32 


12.8 


1*3.6 


1.58 


2-88 


191*5 


20 


21*80 


9 


15/32 


9.81* 


32.3 


1.59 


2-89 


191*5 


20 


21*80 


9 


15/32 


12.1 


1*3.0 


1.1*8 


2-90 


191*5 


20 


21*80 


9 


15/32 


5.U8 


1*0.8 


.9h 


2-91 


191*5 


20 


21*80 


9 


15/32 


6.61* 


25.3 


1.1*1 


2-92 


191*5 


32 


21+80 


9 


15/32 


13.9 


1*1.9 


1.70 


2-93 


191*6 


35 


21*80 


9 


15/32 


9.02 


1*5.1* 


1.17 


2-9U 


191*6 


31* 


21*80 


9 


15/32 


11.0 


1*9.2 


1.26 


2-95 


191*7 


31 


21*80 


9 


15/32 


11*.5 


73.6 


1.16 


2-96 


191*1* 


9 


21*80 


9 


15/32 


5.91 


37.2 


1.02 


2-97 


191*1* 


10 


2I48O 


9 


15/32 


18.1 


1*0.5 


2.33 


2-98 


191*5 


10 


21*80 


9 


15/32 


17.1 


53.3 


1.65 


2-99 


19U5 


10 


21*80 


9 


15/32 


32.6 


61.8 


2.95 


2-100 


19l*5 


10 


21*80 


9 


15/32 


2l*.l 


66.2 


1.87 



404091—56 X 
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Table A-2. Summary ball-bearing data for company B, with computed values for L Wj L 50 , and Weibull slope e — Continued 



Record 


Year 


Number 




z 


D 






e 


No. 


of 


in test 


Load 


Number 


Ball diam. 


^LO 


L *« 


Weibull 




best 


group 




of balls 




5o 


slope 








lb 




in. 








2-101 


19l*5 


20 


21*80 


9 


15/32 


36.1 


71.6 


2.75 


2-102 


19U6 


20 


21*80 


9 


15/32 


63.3 


101* 


3.82 


2-103 


19U6 


12 


21*80 


9 


15/32 


ll*.l* 


59.0 


1.33 


2-101* 


19U6 


11 


21*80 


9 


15/32 


15.1 


92.9 


1.01* 


2-105 


191*5 


10 


21*80 


9 


15/32 


18.8 


39.1* 


2.55 


2-106 


1950 


12 


2U80 


9 


15/32 


5.63 


31*. 7 


i.ol* 


2-107 


1950 


12 


21*80 


9 


15/32 


7.23 


3U.5 


1.21 


2-108 


1951 


30 


21*80 


9 


15/32 


16.7 


71.8 


1.29 


2-109 


1951 


63 


21*80 


9 


15/32 


26.5 


90.3 


1.51* 


2-110 


1950 


23 


21*80 


9 


15/32 


8.35 


1*9.1 


1.06 


2-111 


191*3 


19 


3250 


9 


15/32 


3.79 


9.30 


2.10 


2-112 


1937 


10 


31*70 


10 


7/16 


9.05 


36.6 


1.35 


2-113 


19U1* 


20 


Uooo 


9 


15/32 


2.98 


7.35 


2.08 


2-lli* 


191*3 


19 


2300 


10 


15/32 


22.5 


73.1* 


1.59 


2-11? 


1938 


10 


2730 


10 


15/32 


3.82 


31.7 


.89 


2-116 


191*6 


22 


2660 


10 


17/32 


6.55 


20.8 


1.63 


2-117 


1910* 


20 


2250 


11 


15/32 


17.5 


61*.3 


1.1*5 


2-118 


191*3 


16 


2300 


11 


15/32 


61.7 


152 


2.10 


2-119 


19U5 


1*8 


281*0 


11 


15/32 


18.6 


1*2.7 


2.27 


2-120 


191*7 


28 


281*0 


11 


15/32 


21.6 


66.3 


1.68 


2-121 


191*7 


8 


281*0 


11 


15/32 


11.9 


39.1 


1.59 


2-0.22 


WU8 


8 


281*0 


11 


15/32 


13.9 


50.6 


1.1*6 


2-123 


±9h3 


19 


3200 


11 


15/32 


7.80 


33.1 


1.30 


2-121* 


191*1* 


28 


1*000 


11 


15/32 


3.55 


13.9 


1.38 


2-125 


191*3 


19 


1*000 


11 


15/32 


9.1*0 


23.1* 


2.06 


2-126 


191*7 


23 


6350 


11 


11A6 


1*.76 


22.7 


1.21 


2-127 


191*1* 


20 


12000 


11 


1-1A6 


3.23 


9.86 


1.69 


2-128 


19U1* 


20 


12000 


11 


1-1A6 


2.62 


9.52 


1.1*6 


2-129 


191*1* 


9 


12700 


8 


1-1/2 


7.89 


39.7 


1.17 


2-130 


191*9 


18 


16500 


11 


1-1A6 


1*.93 


20.1* 


1.33 


2-131 


1950 


20 


l65oo 


11 


1-1A6 


6.26 


16.2 


1.98 


2-132 


1938 


8 


565 


7 


5A6 


37.3 


103 


1.85 


2-133 


191*1* 


20 


900 


7 


5/16 


11*.0 


38.6 


1.86 


2-13U 


1938 


10 


i65o 


8 


13/32 


30.3 


87.6 


1.77 


2-135 


191*1* 


20 


2250 


8 


15/32 


25.7 


71.2 


1.85 


2-136 


191*3 


20 


2300 


8 


15/32 


10.5 


60.1* 


1.07 


2-137 


191*1* 


19 


3200 


8 


15/32 


10.3 


2l*.l 


2.21 


2-138 


191*1* 


19 


ItOOO 


8 


15/32 


1*.56 


12.9 


1.81 


2-139 


1937 


10 


1710 


8 


17/32 


25.1 


271* 


.79 


2-11*0 


1938 


9 


2360 


8 


17/32 


1*8.8 


261* 


1.12 


2-11*1 


1937 


10 


2680 


8 


17/32 


7.53 


60.7 


.90 


2-11*2 


1937 


11 


3850 


8 


17/32 


ll*.9 


62.6 


1.32 


2-11*3 


191*7 


21 


7760 


8 


29/32 


1*.57 


U3.1* 


.81* 


2-11*1* 


191*3 


12 


9550 


8 


1-1/16 


3.90 


1*0.7 


.80 


2-11*5 


191*7 


21 


9750 


8 


1-1A6 


15.5 


79.1* 


1.16 


2-11*6 


19U8 


16 


111*00 


8 


1-3A6 


10.2 


U3.9 


1.29 


2-11*7 


191*8 


20 


111*00 


8 


1-3A6 


1*.71 


16.9 


1.1*8 


2-11*8 


191*7 


18 


111*20 


8 


1-3A6 


10.1 


31*. 2 


1.55 
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Table A-3. Summary ball-bearing data for company C, with computed values for Li , L50, and Weibull slope e 



Record 


Year 


Number 




z 


D a 






e 


Ho. 


of 
test 


in test 
group 


Load 


Number 
of balls 


Ball diam. 


L 10 


L 5o 


Weibull 
slope 








lb 




in. 








3- 1 


19h2 


9k 


1580 


7 


9/16 


16.9 


61*. 8 


1.1*0 


3- 2 


19l*9 


29 


790 


7 


9/16 


211 


729 


1.52 


3- 3 


19h9 


35 


1185 


7 


9/16 


7k.h 


287 


1.1*0 


3- h 


19hP 


29 


1600 


9 


1/2 


9.62 


1*0.1 


1.32 


3- 5 


191*5 


10 


1600 


8 


15/32 


11.9 


66.3 


1.10 


3- 6 


191*3 


9 


2275 


7 


17/32 


13.8 


58.0 


1.31 


3- 7 


191*6 


13 


251*0 


8 


15/32 


2.38 


11.3 


1.21 


3- 8 


19W 


12 


25U0 


8 


15/32 


2.38 


11.5 


1.19 


3- 9 


191*9 


12 


1580 


7 


9/16 


8.75 


62.2 


0.96 


3-10 


191*9 


12 


1580 


7 


9/16 


25.7 


113 


1.27 


3-11 


191*7 


21* 


1600 


9 


1/2 


liu5 


113 


0.92 


3-12 


191*9 


12 


610 


8 


5A6 


26.8 


65.6 


2.10 



Table A-4. Summary ball-bearing data for company D, with computed values for Li , L 50 , and Weibull slope e 



Record 


Year 


Humber 




Z 


D a 






e 


No. 


of 


in test 


Load 


Number 


Ball diam. 


L* 
10 


L* 

L 50 


Weibull 




test 


group 




of balls 






slope 








lb 




is. 








li-1 


191*6 


19 


1750 


9 


7/16 


159 


963 


1.05 


Jj-2 


1951 


31* 


1750 


Q 


7/16 


71.7 


526 


0.91* 


k-3 


1951 


56 


1750 


9 


7A6 


113 


582 


1.15 



Life estimates are in hours. 



Table A-5. Summary of test groups of ball-bearing data 



Company 


Number of test 
groups 


Total number of 

bearings in test 

group 


A 


50 
148 

37 
94 
17 
12 
3 


1,259 

3,289 

~~291 
109 


B 


Type B-l 


Type B-2 

Type B-3 _ 


C__ 


D* _ 


Total (all companies) 


213 


4,948 
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*These data were not used in the main analyses 
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SPECIMEN WORKSHEET 



Reference No. 



Bearing Mfg. by- 


Bearing Tested by- 


Date of Test 2-26-46 


Bearing No. 


Load 5CC R.L. 


Speed 2CC0 r.p.m. 


Lubrication: Type Jet Oil 


Frequency- 
Ball No. and Dia. 9 - 1/4" 


Contact Angle U 


Groove Radius: Inner Ring 


51.6* 


Outer Ring 
Number of Rows 1 


55.07S 



Bore 


20 mm. 




O.D. 


42 mm. 




Lot Size 


25 


Taken on 23 



Bearing temperature measured on outer 

ring at point of maximum load 

Material: Type 



Source 

Rockwell Hardness of: 

Inner Ring 63.5 

Outer Ring 64. 

Balls 



Ball Failure 

Inner Ring Failure_ 
Outer Ring Failure^ 



JJL 



J52SL 



2Q% 



Table Ordered According to 
Endurance Life 



JJL 



Brg. 


Endurance 


Type of 


No._ 


Mill. Revs. 


Failure Remarks 


16 


17.88 


Ball 


10 


28.92 


Ball 


5 


33. CO 


Ball 


19 


41.52 


I.R. 


9 


42.12 


Ball 


11 


45.60 


Ball 


15 


48.48 


Ball 


12 


51.84 


Ball 


20 


51.96 


Ball 


18 


54.12 


I.R. 


13 


55.56 


I.R. 


1 


67.80 


Ball 


— 3- 


&rVS9 


— L»itoro Omit tod 


— 3- 

4 


&hQ9 

68.64 


— LjBoro Omit tod 

Ball 


6 


68.64 


L.Bore 


25 


68.88 — > 


Disc. 


22 


84.12 


Ball 


17 


93.12 


Ball 


7 


98.64 


I.R. 


23 


105.12 


I.R. 


24 


105.84 -> 


Disc. 


21 


127.92 


Ball 


8 


128^)4 


O.R. 


14 


173.40 -> 


Disc. 



Test life in 10° 

Median 

Mean 

B-10 

Slope of Curve 

Test No. 

Lot 



revolutions: 



71. 
29. 



2.23 



3183 
71 
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5. Appendix B. Evaluation of L 10 , L 50? and Weibull Slope e, by Using Order 

Statistics for Censored Data 

This is a technical appendix that gives the mathematical basis for estimating, for each 
test group, the values of L w and L 50 for use in the regression analysis discussed in appendix C, 
and also the Weibull slope e. 

5.1. Weibull Distribution 

a. Characteristics 

As noted in the text, the basic assumption for estimating L 1Q , Z 50 , and e for each test group 
was that the probability distribution of fatigue lives of individual bearings could be represented 
by a " Weibull distribution/ ' 8 This means that the observed fatigue lives (number of revolu- 
tions) of all the bearings in a test group of, say, n bearings constitute a random sample of n 
independent observations from a distribution whose cumulative (from above) distribution 
function (hereafter denoted by cdf) is 9 

S(L)=Prob{life>i} 

= exp[-(Z/a) # ], 0<Z<co, (Bl) 

where a and e are the two parameters to be fitted. They are related to L m and L 50 by eq (B2a) 
below. The function S(L) is also termed the "survivorship" function. This distribution is 
one of three limiting types to which the distribution of the smallest member of a sample, under 
general conditions, tends as the sample size is increased indefinitely. (Another type is dis- 
cussed in the following section.) This matter was first studied chiefly by Fisher and Tippett 
[5], and for this reason the type (Bl) is sometimes referred to as Fisher-Tippett type III for 
smallest values. 

There are both theoretical and practical reasons for choosing the Weibull distribution 
(Bl) as the underlying probability distribution for fatigue life. 

Theoretical. Here it is assumed that fatigue is an "extreme-value" phenomenon, related 
in some manner to the strength at the weakest point in the material under stress. The theo- 
retical reasoning that proceeds from this assumption is mentioned by a number of authors, 
and is given explicitly, for example, by Freudenthal and Gumbel in [6, p. 316 to 318]. It 
leads precisely to the form (Bl) (see eq (2.9) in [6]). It is recognized that this statistical 
assumption has not received universal acceptance. This paper is, however, not concerned 
with the relative merits of various statistical theories of fatigue, but merely with consequences 
of a reasonable choice from among them. 

Practical. Application of the Weibull distribution received extensive attention by W. 
Weibull in [19], where he showed that a distribution of the general type (Bl) represented certain 
fatigue-life data quite satisfactorily. In addition, inspection of the special "Weibull" plots 
accompanying the worksheets suggests that many can be fitted satisfactorily by a straight 
line representing a Weibull distribution, as explained below. 

The manner in which these graphs are constructed is described by Weibull in [19]. A 
sample of Weibull-function coordinate paper used for this purpose is included in appendix A. 
The essence of the method is that eq (Bl) may be converted, by taking logarithms twice, into 

e(ln L) - (e In a) =ln[m(l/£)], (B2) 

where "In" denotes the natural logarithm (base denoted by e) and S=S(L). From eq (Bl) 
and the definitions of L w and Z 50 , when L=L i0 , S(L) = . 90; and when L=L 5Q , S(L) = . 50. 
These values substituted in eq (B2) give 

8 So named for W. Weibull (cf. [18], p. 16 fT.), who is considered to be one of the first to study it extensively. 

9 The use of a continuous instead of discrete probability distribution will introduce no appreciable error. 
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e(\n L l0 )-e(h\ a)=ln[ln(l/.90)]= -2.25037 

(B2a) 
Kin L n )—eQn a)=ln[ln(l/.50)]= -0.36651,. 

the values on the right-hand side being obtained from [17, table 2]. These are the relation- 
ships between the parameters a, e, Li , and L 50 . The right-hand numerical values will later 
be denoted by y J0 , y m50 , respectively. Equation (B2) may be written 

ex— a'=y, 
where 

sr=ln L, a'=e In a, y=lnpn(l/S)]. (B3) 

The variables x, y correspond to the two scales shown on the Weibull-f unction coordinate paper 
in appendix A. The variable x, with unrestricted values, corresponds to the horizontal scale 
"Bearing life," having a logarithmic scale. The variable y is represented through the per- 
centage surviving, S, or rather through the (vertical) scale for "bearings tested — percent" = 
percent f ailed 10 =1—S=P, which can vary only between and 1. This scale also has non- 
uniform graduations, given by the iterated logarithm in (B3). 

The Weibull distribution is thus seen to be equivalent to a straight-line relationship, with 
"Weibull slope" e, between the logarithm of fatigue life and an associated quantity y depend- 
ing only on its relative rank when the fatigue lives are arranged in ascending order. Thus, 
goodness of fit of the straight line (B3) is equivalent to goodness of fit of a Weibull distribu- 
tion to the fatigue lives L of an individual test group. In fact, one common method of sta- 
tistical analysis of fatigue-life data (Freudenthal and Gumbel [6]) depends upon the use of 
the classical method of least squares for fitting this straight line. This method is, however, 
subject to certain limitations described below. Instead, an alternative method, presented in 
the following sections, is preferred that fits the distribution of x=h\ L directly by use of order 
statistics. 

b. Limitations of Fitting by Least Squares 

In the classical method of least squares for fitting the straight-line relationship (B3) to 
a test group of ball-bearing data, pairs of values (x t , y t ), i=l, . . . , n, are required. The 
values of x=\n L are obtained from the given data. However, the variable 1 y, measured 
through the percentage failing, P—l—S, presents difficulties. The problem of how to plot 
P is known as the problem of "plotting position." 

It seems clear that the values, P iy of the plotted variable, P, must somehow be related to 
the rank order of the bearings as they fail. A natural choice is the percentage failing: P=f/n, 
where / is the rank order of failure in a test group of n. This is not advisable for reasons dis- 
cussed at length by Gumbel in [9, p. 14], where he advocates the plotting position //(/i+l). 11 
Other workers take different positions, and the question of plotting position must be regarded 
as still unsettled. 

A second difficulty with the use of least squares is that as usually used it fails to take 
adequate account of the number of items remaining intact ("runouts") in the incompleted 
tests. As a final point, it is to be noted that the successive plotted points are not independent, 
as they represent the observed lives in increasing order. A correct use of least squares pro- 
cedures would have to take into account all the intercorrelations, which is not done in the 
usual application of the "method of least squares. " The method of order statistics described 
in section 3 has the advantage of avoiding the above limitations of the least squares method. 



10 The symbol P as used here should not be confused with the same symbol for load used in the life formula. In any event, the meaning 
will be clear from the context. 

11 This plotting position was also used by Weibull in [18] (cf. eq (72) and the vertical scale in figures 3 and 4 therein). 
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5.2. The Extreme- Value Distribution 

a. Relation to Wei bull Distribution 

The preceding section indicates that logarithms of lives, rather than lives themselves, are 
the natural units in which to carry out the analysis. This idea has also been adopted by those 
who do not use the Weibull distribution, either because they are unaware of its existence or 
because they do not feel it fits their data. 

If the Weibull distribution is adopted for fatigue life, L, then the variate, z=ln L, has 
the nonnormal cumulative distribution function 12 

0(x) = Prob { In (life) > x } = Prob { life > i } 

= S(e x ) = exp(-a- e e ex )=exp[- e (x - lna) /e ~ 1 ], - oo <>< o> . 

This may be written 

G(z)=*(y) = exp(-*), (B5) 

where 

y=(x—u)!l3, — ocx^aK^ 00 , (B6) 

and 

tt=]n a, 0=1/6 (B7) 

are the two parameters. The distribution, <£(?/), considered as a distribution of the "reduced 
variable/' y, has standardized parameters u=0, 0=1, and is called the "reduced distribution." 

The form (B5) is another of the three asymptotic distributions of extreme values, some- 
times designated as Fisher-Tippett type I for smallest values. This distribution has been 
studied extensively., chiefly by E. J. Gumbel (e. g., [7, 8, 9]). In this paper, the term "extreme- 
value distribution" will be given to the distribution of smallest values (unless otherwise speci- 
fied), although this name is frequently given to the largest-values case. 

From the above discussion, it is apparent that methods pertinent to the type I extreme- 
value distribution (B5) are appropriate. For this purpose there is available a mathematical 
approach recently developed by one of the authors of this report, and described in detail in [13]. 

b. Characteristics 

A description of the extreme-value distribution (B5), together with an interpretation of its 
parameters in terms of life estimates (or rather their logarithms), is essential to an understand- 
ing of the application of the method of order statistics in this paper. It will be seen that the 
problem of estimating life is equivalent to that of estimating the parameters u and 0. 

The parameters of the extreme-value distribution (B5) are depicted in figure 1 (page 291). 
The quantity u is the position of the mode or highest point of the (frequency) distribution. 
The quantity is a scale parameter, analogous to the standard deviation, a, in the case of the 
normal distribution. In fact, is t/q/tt (about %) times the standard deviation of the extreme- 
value distribution. 

Although the two parameters, u, 0, completely specify the distribution, it is very useful 
to introduce related quantities of the form 

t=u+py, (B8) 

which are linear combinations of parameters u and and may thus also be regarded as param- 
eters when known values are later assigned to y. Introduction of t makes it possible to esti- 
mate u and simultaneously. Thus if t can be obtained as a-{-by with o and b known and y 
arbitrary, then the values u=a, 0=6 can be read off at once. 

The parameter t has another highly important meaning. In figure 1 the area F under the 
distribution to the right of the ordinate erected at t represents the probability that a value 

12 Cf. Freudenthal and Gumbel [6], eq (2.5), (2.6), (2.8), (2.9). 
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X IO = lnL IO =U ^^y.90 



yqo=- 2 - 25037 



x 50 =lnL 50 =u+ i9y 




y. = -0.36651 



t F = ln L 50 - 



x = lnL 

Figure 1. General form of extreme-value distribution (for smallest values) showing relationship of 
parameters tF, Xio = ln Li , and Xw = ln L 50 , to u and 8. 
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larger than t will occur. Thus t is a function of F and may be written t F , as shown ; it is desig- 
nated the "upper lOOi^-percentage point" of the distribution. For example, if F=.90, then 
t=t,go represents a value of x=\n L, which will be exceeded by 90 percent of the population. 
This is associated with rating life L 10 (life exceeded by 90 percent of bearings) by the relation 

£.9o=£io=ln L i0 , (B9) 

where x represents life in logarithmic units. Similarly, for median life, 

£.5o=%)=ln L 50 . (BIO) 

Since the t's are regarded as parameters of the distribution, so also are x 10 and x 50} and therefore 
Zio and L 50 . These are not, of course, all independent. 

In general, we have the percentage point t F , which, expressed in terms of the original param- 
eters u and ft, may be written in the form (B8) : 

t F =u+/3y F} (B8a) 

where y is a quantity depending only on the probability F, determined as follows. We have 
from (B8a) 

y F =(tr-u)l0, (Bll) 

i. e., y F is the value of (x—u)/($ when x takes the value t F . But by definition of the probability F, 
in view of (B5), (B6), and (Bll), 

F=Froh{x>t F } = G(t F ) = <f>(y F ) = exp(-e VF ). (B12) 

Thus, solving for y F , we obtain 

y F =ln(-lnF). (B13) 

This is the reduced variable corresponding to the probability F, and may be obtained by a simple 
change in sign from table 2 of [17], which tabulates the function 

-ln(— In %), 

where $„, a probability, takes on values from to 1. Thus, 

for F=.90, y F = -2.25037; 



for F=.50, y F = -0.36651. 



} (B14) 



The above discussion shows that both x 10 and x 50 (rating and median lives in logarithmic 
units) may be determined once the general percentage point (B8a) is estimated by giving the 
two particular values (B14) to y F . 

c. Conversion From Largest to Smallest Values 

The methods and numerical results developed in [13] were for problems, such as maximum 
gust-loads on airplanes, that required the distribution of largest sample values. In order to 
adapt this material to the distribution of smallest values (B5) required here, the relationships 
of symmetry involved in the reversal of direction must be examined with care. To avoid 
confusion, it is necessary to use subscripts L and S to distinguish between quantities related 
to the largest-value distribution from those related to the smallest-values case. No generality 
is lost by use of reduced variates. Thus, in (B5), x will be replaced by the reduced variate 
y, and, for simplicity, the symbol 0(y) will be used instead of $(y) : 

G(y) = $(y) = exp{-e«). (BIS) 
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From this, the ("cumulative from above") distribution of smallest values is 

Froh{Y s >y}=G s (y) = exp(-e«), -«<y<«>, (B16) 

where Y s denotes the reduced smallest value. The corresponding distribution of largest 
values is (see Gumbel [9, eq (I), p. 21]) 

Frob{Y L >y}^H L (y) = l-Frob{Y L <y} = l-exp(-e^) = l-Gs(-y) f -co<</<oo, (B17) 

from (B16). 

The corresponding relation for the density functions is obtained by differentiation, with 
gs(y) = Gs(y), and h L (y) =H' L (y) : 

9s(y)=h L (-y). (BIS) 

Hence the two distributions are merely mirror images of each other. The moments of the 
distributions are related as follows : 

v ks ^E s (f)=£ gi y k gs(y)dy=j^ {-yJh L {y f )dy=(-lfv kL . (B19) 

Thus, the means differ in sign and the variances are identical: 

vm= -y=-viL, (B20) 



*%=^=*b (B21) 

These values are given, for example, in [9, p. 23, eq (3.27)]. 

Finally, we need the relationships between moments of the order statistics for the two 
distributions. As the smallest-value distribution is a reversal of the largest-value distribution, 
it is natural to reverse the arrangement of the order statistics as well. This gives simpler 
results. Thus we are interested in the ith order statistic in the series 

OS): y[>y' 2 > . . . >y' t > . . . >y' n , (B22) 

where the parent distribution is that of smallest values. Primes will be used as a reminder 
that the order is descending, not ascending. Thus in tables B-2 and B-3 the absence of primes 
indicates that the order statistics are in increasing order. 
(B22) is the analogue of the series 

(£): Vi<V*< . . . <yt< ...<». (B23) 

of order statistics for the largest-value parent distribution. Whenever a distinction is neces- 
sary the subscripts S or L will be used with the y's. 

From (B18) it may seem intuitively (and may be justified rigorously) that the distribu- 
tions and moments of the order statistics follow the same symmetry relationships as the parent 
distributions, namely, 

E s (y'r)=(-l) m E L (yT) 



(Ts,ij—< T s(yi,y'j)=<rL(yi,yj)=(rL,ij- > 



(B24) 



In other words, the even moments remain the same; the odd moments change only in sign. 

The above development shows that the numerical results for moments of order statistics 
previously obtained in [13] for the largest-value case can be used here for smallest values 
without any substantive change. 
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5.3. Method of Order Statistics for Censored Samples 

a. For Small Samples 

Consider an independent random sample of n items from the distribution of smallest 
values, of which only the k smallest values can be observed. In view of the preceding dis- 
cussion, it is desirable in the theoretical development to deal with the order statistics in descend- 
ing order: 

(x[>X 2 > . . . >X n _ k )>X n _ k + y> . . . >Zn, (B25) 

where the parentheses denote the (n—k) (largest) unobservable values, and the remaining k 
values are known. This arrangement materially simplifies the exposition. Primes will again 
be used to denote descending order to distinguish from ascending order, which will occur in the 
later parts of this section. 

From the k known values it is desired to determine an estimator 

Tn, k -=w[x f n . k+l +W2Xn^ k+ 2+ . . . +w' k x' n , k<n, (B26) 

(i. e., the weights w'j) of the general parameter, 

t F =u+fiy F , (B27) 

of the extreme-value population (B5), such that T' in (B26) is (1) unbiased and (2) of minimum 
variance. Mathematically, this means that 

E(T') = t F , (B28) 

where E denotes mathematical expectation, and 

Var (T') = a minimum, (B29) 

subject to the above condition. 
From (B6), 

x=u+Py, (B30) 

where y is the reduced variable and x the observed variable. From this the following relations 
for the order statistics x t and y t are apparent: 

Xj=u-\-$y'u j=n—k+lj n—k+2, . . . } n f (B31) 

Xn-Jc+l >Xn-Jc+2> • • -^!X n , (B32) 

^n-*+l>2/a-*+2>. • ->Vn, (B33) 

Esiz'^u+pEstf,). (B34) 

The values E s (y'j) may be obtained with the aid of the table in [14]. This table gives the 
values of E L (y' r ) where the order statistics, y' T , are in descending order (as indicated by the 
prime). The means needed in (B34) are obtained from (B24) and the evident symmetry 



relations 

E L {yf) = {-\rE L {y' n m -i + x) 




as 




E s (y' t ) = -E L (y' H „ s+1 ). 


(B35) 


Reference [14] gives the values of E L (y' r ) for r=l(l)min(n,26), « = 1(1)10(5)60(10)100. 




From (B28) and (B31), 




E(T')=j:w'i [u+pE(y' n _ k+J )] = t F =u +0y F . 

7 = 1 


(B36) 
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This is required to be an identity for all values of the parameters u, fi. Equating their co- 
efficients gives the two conditions on the weights, w'ji 

k k 

I>;=1, S[tf(yi-* + ,)]ti>;=!fr, (B37) 

j=l j=l 

where the numerical values E(y' n _ k+j ) may be obtained from [14] as already indicated. 
For the variance condition (B29), we have, in view of (B26), 

k k k 

^i(2 7, )=S^ V(^_ fc+ ^)+S S^i^i^-fc+i, x'n-k+j)- (B38) 

VaKrO^lS^Vl^.+ ZSV^^;.^,,^^^ 2 -^^'^ (B39) 

= minimum subject to (B28). 

Use of Lagrange multipliers in the same manner as in [13, pp. 50-52] gives, after differentiation, 
the conditions on the weights: 






(B40) 



For each fixed value of k<n there arc k linear equations which, with the two in (B37), form a 
simultaneous system of (£ + 2) equations in the (k-\-2) unknowns, w[, w' 2 , . . ., w' k , X, /x. The 
values of X and /x are useful as a check, because, if (B40) is multiplied by Wj and summed for 
i=l,2,. . ., k, the result is, in view of (B37) and (B39), 



Vl% in +\+nY F =0. 



(B41) 



The minimum value, V*% n , will be denoted by Qn.*« 

In general, there will be a set of (i+2) linear equations to solve for each k = 2, . . ., n. 
(1) Case k=n. For k=n, the matrix of coefficients and right-hand "constant terms" of 

(B40) and (B37) is the (w+2) by + 3) matrix 

1 Ey[ 
1 Ey' 2 



A° = 



022 



°"w2 



0"2ra 



#2/» 



1 


1 


1 








l 


Ey[ 


^2 • 


Ey'n 








Vf 



(B42) 



The ordinary (n + 2) by (n + 2) matrix of coefficients, without the constant terms, will be de- 
noted by A n . If T n denotes the vector column of constant terms, then 

Aj=[Adrj, (B43) 

and the linear system of (ft + 2) equations may be denoted 

A n W:=T n , (B44) 

where W' n denotes the column vector of the (n+2) unknowns w[, w' 2 , . . ., w' n , X, \x. 
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The coefficients of the unknowns in (B44) involve the means E(y'j), already discussed, and 
covariances a[j. These values are given in table B-l for n=2 to 6. The a'a were computed 
by the method developed in [12]. Table B-l also indicates how the moments for the largest 
values case can be obtained simply from those shown. 

The (n-\-2) solutions of (B44) are all expressible linearly in terms of the components of T n . 
Thus the solutions all take the form 



wj=aj+%jr, i=l,2, 
\=c[+d[y F 

fjL=c' 2 +d2y F 



(B45) 



Substituting these w) in (B39) gives an expression of the form 

Qn,n=Y^=A f n +2B , n y F +C , n yl. (B46) 

The quantities a' h V 6 for the weights w' h for n = 2 to 6, are shown in table B-2. The coefficients 
A' n , B' n , C' n of Qn, n, and the values of Q' ntU evaluated at 2^= .90, .50, corresponding to L 10j L 50 . 
respectively, are given in table B-3. 

Calculations were limited to n = Q in this paper, in view of the diminishing returns in 
"efficiency" (see below) for increasingly larger amounts of computing. Methods suitable for 
larger values of n are discussed later. 

Table B-3 shows that as sample size increases from n=2 to 6 (in the case k=ri), the 
variance diminishes for the percentage-point parameters t F for F=.90 and .50, i. e., t F =Xi =\n 
Zio and t F =x 50 =\n L 50 . This is a common characteristic of the behavior of estimators for 
increasing sample size. Another method whereby estimators may be compared is through 
their efficiency. 

Efficiency is a measure intended to provide a convenient standard of comparison for 
estimators. This is done for two estimators to be compared by dividing the variance of each 
into a theoretical "smallest" variance, Q LB , known as the "Cramer-Rao lower bound." Further 
details in the case of complete samples where k=n, as here, may be found in [13, p. 14 and 
15]; values of Q LB are also indicated in this reference in table III (a). 

Table B-4 shows the efficiency values so obtained, for the case k—n, n = 2 to 6, as regards 
the order-statistics estimators for the parameters Xi =ln £ 10 , x 50 =ln L m . 

These values show that for x i0 , the efficiency, starting with under 70 percent for n=2, 
increases rapidly until 89 percent, out of a possible maximum 100 percent, is reached for n=6. 
A 90-percent-efficient estimator is generally considered to be quite good. As regards x 50 , 
the efficiency is well above 95 percent for all the values of n, and for n = 6 exceeds 99 percent. 
In view of results of this nature, and because of the increasingly heavy computations necessary, 
calculations were not carried beyond n=§ in [13]. 

The above applies to estimation of the parameters x i0 and x 50 , which it will be recalled are 
the logarithms of the actual life estimates Z 10 , Z 50 . It is believed that efficiency of the method 
of order statistics in obtaining estimates of actual life Z 10 , Z 50 is probably reasonable, in view 
of its high efficiency in estimating the logarithms, #io, x 50 . 

(2) Case k<^n. For the case k<^n, the procedure is very similar. One starts with a 
(&+2) by (& + 3) order matrix A? derived from A° in (B42) by striking out the first (n—k) 
rows and columns. One proceeds in this manner for k=n—l, n—2, etc., until when k=2 
the matrix becomes 



AS= 



&n-l,n-l 



G n, n — \ 



Vn-\,n 



1 




Jty'n-i Ey'n 

representing a set of 4 equations in 4 unknowns. 



Ey'n-i 
Ey'n 







1 
Vf. 



(B47) 
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The resulting weights w) and variances Q r Ufk were obtained in similar [manner to those 
for k=n in (B45) and (B46). These, it will be recalled, are primed quantities, associated with 
descending order of the order statistics. Because the observations, x, for successive failures 
naturally occur in ascending order, it is more useful for actual application, in contrast to theo- 
retical development, to tabulate the weights and covariances for the order statistics in ascending 
order. This has been done in table B-2, giving the weights Wi=a i J rbiy F , and in table B-3, 
giving the variances Q ntk =A+2By F J rCy F for the estimators T n>k formed with the above 
weights. These variances are also evaluated for the parameters #i =ln i 10 , ^ 50 =ln L 50 . The 
relationships of these unprimed quantities to the primed ones of the previous theoretical 
development is merely a reversal of the order throughout, as indicated by subscripts: i. e., every 
a'i is changed to the corresponding a k _ i j rX and similarly for b \ and w\. The variances Q may 
be shown to remain unchanged. 

b. Extension to Larger Samples 

Samples of more than six items are broken up into independent samples of 6 with a re- 
mainder subgroup, if necessary, of from 2 to 5 items. Because the endurance data were arranged 
in increasing order of life, independent random subgroups could not be obtained by simply 
taking groups of 6 in the (numerical) order in which they appeared on the worksheets. It was 
therefore first necessary to randomize the endurance lives on each data worksheet. This was 
accomplished by use of random numbers that were generated in the electronic computer (the 
SEAC) as needed. 

Such artificial randomization is not desirable when it can be avoided, because the results 
of the calculations are then not unique, but may depend to a limited degree on the particular 
set of random numbers used. 13 It is therefore recommended that when the bearings in a test 
group are to be simultaneously run on a battery of fatigue-testing machines, the individual 
bearings should be recorded in advance in some more or less natural order independent of the 
order in which failure takes place in the course of the test. Natural order might be order of 
manufacture, order of testing, etc. 

In the present investigation, each subgroup was treated as a random sample by the methods 
already developed for size 6 or less. That is, a "subestimator" was calculated for each sub- 
group and the results averaged to produce an over-all sample estimator. 

An estimator, both for the individual subgroup and for the over-all sample, was obtained 
for each of the four population quantities: 

u, j8, t. QQ =x 10 =ln L l0 = / u+y. 9O p j t M =x 50 =h\ L 50 =u+y M P- (B48) 

For subgroups, these four parameter estimates are given by (caret denotes "estimate of") 

A k A k A A A A A 

u^^diXt, |3=SMo x 10 =T ntk (10)=u+y^, x 5 o=T nik (50)=u+y m50 p y (B49) 

t=l i=l 

where //„„ —2.25037, y.m= — 0.36651, and&i<a? 2 ^ . . . <x k , 2<k<n<6, are the logarithms 
of the actual observed lives in a subgroup arranged in ascending order, and the a t and 6 7 are 
read directly from table B-2. For the over-all sample estimator, the subestimators T n>k are 
merely averaged. 

For later use (appendix C) the variance of the over-all estimator, T, and its relation to 
sample size will be considered here. Consider first the case of a complete sample, where no 
intact bearings arc present because the test is run to completion. Let n be the sample size; 
then there are two cases, according as (1) n<6, or (2) n>6. 



3 This effect can be reduced somewhat by making a duplicate run and averaging the results, as was done in this study. 
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Figure 2. Relationship of variances Q lQ and Q 5Q to sample 
size n for n=2 to 6 (logarithmic scale in each direction). 

Qio is variance of estimator of £io=ln Zm 

Q50 is variance of estimator of Jso=ln Lao 

All Q's are in units of 2 

(1) n<6. Table B-3 gives the numerical variances, Qi and Q 50 , for n = 2 to 6. These 
values are plotted in figure 2 on double-logarithmic paper. The values for Q 50 (right-hand 
scale) are seen to lie on a straight line of slope negative unity. This shows that at least in this 
case, variance is inversely proportional to sample size. For the other case, Qi , a straight line 
also gives a reasonably good fit, and its slope appears to differ only a little from — 1 . Hence 
the underlined statement is approximately applicable here too. 

(2) n^>6. If a sample of size n^>6 is broken into equal subgroups (of size 6, for example) 
there will generally be a remainder of size less than 6. The preceding development, when 
suitably modified, shows that for large n the influence of this remainder is small compared to 
the remaining bulk of the sample and thus the rule in question holds approximately in this case. 
Agreement with the rule is less close for a few cases of moderate n, but for simplicity the inverse 
relationship will be taken as a reasonable rule of thumb in all cases for the over-all purposes 
of analysis. 

Two complete runs were made on the SEAC for each of the 213 test groups of data, and 
the two results were averaged for each group, giving values of the averages 



u, ft ?.9o=ln iio, ^.5o=ln L 5{ 



(B50) 



From these, the values of L m and L 50 were obtained from a table of exponentials and the 
Weibull slope e=l/p obtained as a consequence of formula (B5). An example showing the 
steps in calculation of L 10 , i 50 , and e is discussed below. 

Because this investigation represents probably the most extensive mass fitting of the 
Weibull distribution made to date, a tabulation of the 213 values of the parameter e will be of 
considerable value to future applications of this distribution. This is shown in table B-6. 
The corresponding histogram is given in figure 3. Particular items of interest are 

mean of e (all 4 companies) = 1.51, 

median e (all 4 companies) = 1.43. 

Note that 50 percent of the values are in the interval 6=1.17 to 1.74. 
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Figure 3. Histogram for the 218 estimates of parameter e for companies A, B, C, and D combined. 

I )ata in table B-fi 

c. Worked Example 

The example that will be given to illustrate the foregoing procedures will be the one that 
was worked, out as a "test problem" for the SEAC before using the full set of data. The test 
group of bearings selected for this purpose was No. 1-1 in table A-l, for company A. The test 
group consisted of 24 ball bearings, of which 4 remained intact when the test was discontinued. 
The details of the computation for obtaining values of Li , L 50 , and Weibull slope e from the 
test group of data are contained in table B-5 and described in the steps below. 

The endurance lives in observed (increasing) order are listed in column (1). The arrows 
indicate the four "run-outs/' or "intacts," whose testing was discontinued at the number of 
million revolutions indicated. All that is known about these four bearings is that their fatigue 
lives exceeded the values shown. 

Step 1. Randomization. The order of endurance lives in column (1) was randomized by 
use of a set of random numbers generated in the SEAC as part of the computation work. The 
result is shown in column (2) of table B-5. 

Step 2. Subgroups. The lives in randomized order were divided, as shown by the lines of 
separation, into subgroups of size n = 6, the maximum size for which the order-statistics weights 
had been computed. 14 Each subgroup was then prepared for the application of the order- 
statistics method by rearranging in increasing order (column (3)). Natural logarithms were 
then taken as in column (4) . 

Step 3. Weights. Each subgroup was regarded as consisting of k actual observations out 
of a censored sample of n. It happened here that n was 6 for every subgroup ; k took the values 
6, 5, 6, 3. These values are shown in the subscripts of T n , k written in column (3), and they 
determined the weights a t and b t to be selected from table B-2. These weights are represented 
in columns (5) and (6) . 

Step 4- Cross-products. The cross-products 






^LjbiX t 



were then evaluated and placed as shown for each of the subgroups. 



14 As sample size 24 is an exact multiple of n=G, it so happened that there was no "remainder subgroup" in this case. This will not usually 
be true, but the procedure is identical for other values of n, differing merely in the numerical weights to be used. 
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Step 5. Estimates. A simple arithmetic average of the four values was taken for each of 
the two columns (5) and (6), and denoted by Ti=u, T 2 =(3, respectively. These are the order- 
statistics estimates of the two parameters u and f$ of the extreme-value distribution that 
represents the underlying Weibull distribution. 

The reciprocal of £ yields the Weibull slope e= 1.32497. (A) 

The following logarithmic life estimates were given by following linear combinations of 

A A 

the estimates u and ft, using the given values of y Mj y, 50 : 

&io=ln L 10 =u—2. 25037 /3=2.982305 
x 5 o=lnZ 5 o=^-0.36651 /3-4.404120 
rating life Zio^antilog (x 10 ) (base e) = 19.2333 million rev (B) 

median life Z 5 o=antilog (x 50 ) (base c) =81.7872 million rev. (C) 

These three values (A), (B), and (C) represent the outcome of the calculation. 

In the full-scale computing program, calculations were carried out by the SEAC to a larger 
number of places than is shown in the table for presentation purposes. In general, however, 
the number of places shown here should be adequate. The values L 10 , L 50 , and e shown here 
differ slightly from those recorded in table A-l because the latter represent averages of two 
separate runs. 

Table B-l. Means, variances, and covariances of order statistics yi in samples of n from the reduced extreme- 
value distribution G(y)=exp( — e~ y ), n = 2 to 6 

For distribution of largest values, y\<y%< . . . <y n 
For distribution of smallest values, y[>y' 2 > . . . >y' n 



n 


i 


Means* 

= -E L ( yi ) 


Variances and covariances,* 0-^=0-^= <r ij = o- ji 


j=l 


i=2 


i=3 


i = 4 


j=5 


7 = 6 


2 


i \ 


0. 11593 152 
- 1. 27036 285 


0. 68402 804 


0. 48045 301 

1. 64493 407 










3 


! J 


0. 40361 359 

-. 45943 263 

- 1. 67582 795 


. 44849 796 


0. 30137 144 
. 65852 235 


0. 24375 810 
. 54629 438 

1. 64493 407 








4 


ii 


0. 57351 263 

-. 10608 352 

-. 81278 175 

-1. 96351 003 


. 34402 417 


. 22455 344 
. 41553 113 


0. 17903 454 
. 33720 966 
. 65180 236 


0. 15388 918 
. 29271 188 
. 57432 356 

1. 64493 407 






5 


( 1 

2 

{ 3 
4 

I 5 


0. 69016 715 

. 10689 454 

-. 42555 061 

- 1. 07093 582 

-2. 18665 358 


. 28486 447 


. 18202 536 
. 30849 748 


. 14358 737 
. 24676 731 
. 40598 292 


0. 12257 865 
.21226 644 
. 35267 072 
. 64907 319 


0. 10901 329 
. 18967 383 
. 31716 095 
. 58991 519 

1. 64493 407 




6 


(! 


0. 77729 368 

. 25453 448 

-. 18838 534 

-. 66271 588 

- 1. 27504 579 

-2. 36897 513 


. 24658 20 


. 15496 74 
. 24854 56 


. 12121 61 
. 19670 62 
. 29761 59 


. 10291 64 
. 16806 28 
. 25616 60 
. 40185 52 


0. 09116 19 
. 14945 32 
. 22887 90 
. 36145 55 
. 64769 96 


0. 08285 42 
. 13619 10 
. 20925 46 
. 33204 51 
. 59985 67 

1. 64493 41 



*The means are for smallest values (denoted by subscript S) ; for largest values, change all signs and reverse 
order of y's. The an are the same for both. 
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Table B-2. Weights Wifor the order-statistics estimator T n ,kfor the "parameter tF = u-\- Puf of the extreme-value 
distribution {smallest values) from a censored sample of n — 2 to 6, where only the k smallest values are known 

T n ,k=wiXi-\-W2X2+ . . . -\-WkXk,Wj=aj+bjy F 

£i<T2< . . . <Xk,k=2ton 



n 


A; 




Xi 


x 2 


Xz 


Xi 


x- 


x« 


2 


2 


ft 


0. 0836269 


0. 9163731 










-. 7213475 


. 7213475 












f 2 


{&*• 


-. 3777001 


1. 3777001 










3 


-. 8221012 


0. 8221012 










[3 


{hi 


. 0879664 


. 2557135 


0. 6563201 










-. 3747251 


-. 2558160 


. 6305411 










( 2 


ft 


-. 7063194 


1. 706319 1 












-. 8690149 


0. 8690149 










4 


3 


ft 


-. 0801057 


. 0604310 


1. 0196741 








-. 4143997 


-. 3258576 


0. 7402573 










L 


{&.• 


. 0713800 


. 1536799 


. 2639426 


0. 5109975 








-. 2487965 


-. 2239192 


-.0859035 


. 5586192 








' 2 


[b] 


. 9598627 


1. 9598627 












. N962840 


0. 8962840 












3 


ft 


-. 2101141 


-. 0860231 


1. 2961372 








5 


-. 4343419 


. 3642463 


0. 7985882 








< 
4 


ft 


-. 0153832 


. 0519642 


. 1520750 


. 8113440 








-. 2730342 


-. 2499 429 


-. 1491094 


. 6720865 








> 5 


ft 


. 0583502 


. 1088236 


. 1676091 


. 246283 1 


0. 4 189341 






-. 1844826 


-. 1816564 


-. 1304534 


-. 0065354 


. 5031278 






[ 2 


ft 


1. 1655650 


2. 1655650 












0. 9141358 


0. 9141358 












3 


ft 


-. 3153968 


-. 203131.") 


1. 5188283 










. 4466018 


-. 3880492 


0. 8352510 








6 


< 4 


ft 


. 0865378 


-. 0280534 


. 0649390 


1. 0496521 






-. 2858647 


-. 2654 739 


-. 1858756 


0. 7372142 








5 


6; 


. 0057311 


. 0465729 


. 1002523 


. 1722784 


. 6751653 






-. 2015431 


-. 1972753 


-. 1536040 


-. 0645894 


. 6170118 






, 6 


[b] 


. 0488669 


. 0835221 


. 1210527 


. 1656192 


. 2254909 


0. 3554481 




-. 1458072 


-. 1495332 


-. 1267277 


-. 0731937 


. 0359868 


.4592751 
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Table B-3. Variance Q n ,kP 2 of order-statistics estimator T Hi k, given in table B-2 1 and its numerical values 
Q n ,k(10)=Qio, Q n ,h(b~0)=Q5o for estimators of parameters t.go=xio = lnLio, t.5o=xso=lnL 50l respectively, for a 
censored sample of n=2 to 6 

Variances in units of /3 2 

Q n , k =A+2By F +Cy F 2* 

x\<X2< . . . <Xk, k=2ton 



n 


k 


A 


B 


C 


o 10 


Qso 


2 


2 


0. 6595467 


0. 0643216 


0. 7118574 


3. 975015 


0. 708021 


3 


{j 


. 9160386 
. 4028637 


. 4682465 
-. 0247719 


. 8183654 
. 3447117 


2. 952920 
2. 260033 


. 682735 
. 467327 


4 


i 1 

1 4 


1. 3340189 
0. 4331573 

. 2934587 


. 7720298 

. 1180273 

-. 0346903 


. 8670220 
. 3922328 

. 2252828 


2. 250056 
1. 888278 
1. 590460 


. 884572 
. 399329 
. 349150 


5 


1! 


1. 7891720 

0. 5293953 

. 2918142 

. 2313953 


1. 0115594 

0. 2353740 

. 0385708 

-. 0339905 


. 8950462 
. 4168155 
. 2537913 
. 1666472 


1. 769068 
1. 580861 
1. 403458 
1. 228307 


1. 167910 

0. 412852 

. 297633 

. 278697 


6 


(1 

I 6 


2. 2440055 

0. 6529409 

. 3237185 

. 2236063 

. 1911738 


1. 2082248 

0. 3332488 

. 1020223 

. 0105329 

-. 0313731 


. 9132926 
. 4321160 
. 2697162 
. 1861069 
. 1319601 


1. 431164 
1. 341381 
1. 230430 
1. 118677 
1. 000644 


1. 481035 

0. 466709 

. 285165 

. 240885 

. 231897 


*For 


Qw,v F 


= -2.25037; for M> 


y, = -0.36651. 









Table B-4. Efficiency of order-statistics es- 
timator of logarithmic life X\o = ln Lio and 
Xoo = In L 50 for complete samples (k = n) of 
size n = 2 to 6 



n 


Efficiency (in percent) 
with respect to 


3io = ln Lio 


x 5 o = \n L 50 


2 
3 
4 
5 
6 


67.2 
78.8 
84.0 
87.0 
89.0 


97. 4 

98. 3 
98. 7 

98. 9 

99. 1 
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Table B-5. Example showing use of order statistics method of computing values of Lio, £50, and e from endurance 

data 



Test group 1-1, company A 



1 


2 


3 


4 


5 


6 


Endurance (million revolutions) 


Natural loga- 
rithms Xi 


Weights 


Observed 
order 


Randomized 
order 


Ascending 

order within 

subgroups 


a { 


bi 


6.0 


119.0 


27.5 


Z!=3. 31419 


0. 048867 


-0. 145807 


8. 6 


138. 


69. 


0-2 = 4. 23411 


. 083522 


-. 149533 


17. 8 


146.0 


119. 


0:3 = 4. 77912 


. 121053 


-. 126728 


18. 


151.0 


138.0 


z 4 = 4. 92725 


. 165619 


-. 073194 


27. 5 


27.5 


146. 


o: 5 = 4. 98361 


. 225491 


. 035987 


33. 5 


69.0 


151. 


OT 6 =5. 01728 


. 355448 


. 459275 






T 




1X0:;= 4. 817310 
1 

0. 005731 


Zte=0. 400992 
1 

-0. 201543 


50.5 


(150. 0)-> 


6. 


o?i=l. 79176 


51. 5 


8.6 


8. 6 


o 2 = 2. 15176 


. 046573 


— . 197275 


69.0 


51.5 


51. 5 


X3 = 3. 94158 


. 100252 


-. 153604 


74. 


89.0 


89.0 


Xi = 4. 48864 


. 172278 


-. 064589 


74. 


109. 


109. 


0:5 = 4. 61935 


. 675165 


. 617012 


89. 


6. 


(150. 0)-* 


(xt) 










T 65 




£a#»=4. 446363 

1 


2>,r/=l. 213655 

1 


109. 


74.0 


18.0 


Zi = 2. 89037 


118.0 


181. 


33. 5 


X2 = 3. 51155 






119. 
138. 


141. 
18.0 


74. 
118. 


o: 3 =4.30407 
0:4 = 4. 77068 


T Weights are same as"| 
L for first subgroup J 


141.0 


33.5 


141. 


0:5 = 4. 94876 






144.0 


144.0 


144. 


o: 6 = 4. 96981 










17. 8 




2>iZ; = 4. 628081 
1 

-0. 315397 


EM»=0. 619440 
1 

-0. 446602 


146. 


17.8 


Oi = 2. 87920 


(150. 0)-> 


(153.0)-* 


50. 5 


0:2 = 3. 92197 


-. 203432 


-. 388649 


151. 


(153. 0)-> 


74. 


o 3 = 4. 30407 


1. 518828 


. 835251 


(153. 0)-* 
(153.0)-* 
(153.0)-* 


(153.0)-* 
50. 5 


(153.0)-* 

(153.0)-* 
(153.0)-* 

T$ § 3 












74. 








IX-o:; = 4. 831197 

1 


£6^=0. 784853 

1 


Over-all e 


stimator _ 






?!=4. 680738 


f 2 = 0. 754735 














SUMMARY 






f\ = 


rt =4. 680738 


T 2 =j8=0. 7547 


35 c=l/j5= 1.32497 






Vm=~ 


-2.25037 y.Bo= 


= -0.36651 






&io=ln Lio= 


A+ 2/ .90/8 = 2.98230 


5, L 10 = 19.7333 






& £ o=ln £50 = 


A-f 2/ .50^=4.40412 


0, Lfio = 81.7872 
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Table B-6. Tabulation of estimated values of Weibull slopes, e, for the 218 test groups of companies A, B, C, and D 



Weibull slope, e 


Company 


Total 


A 


B 


C 


D 


0.50 to 0.74 


2 

6 

16 
11 

7 
6 

2 


1 

11 

22 

40 
32 
17 

12 
2 
4 
4 

1 






3 

20 

43 
56 
40 
23 

15 
2 
4 
4 

1 


.75 to .99 

1.00 to 1.24__ 


2 

3 
5 

1 


1 
2 


1.25 to 1.49 


1.50 to 1.74 




1.75 to 1.99 




2.00 to 2.24 ___ 


1 




2.25 to 2.49 

2.50 to 2.74 

2.75 to 2.99 












3.00 to 3.24 __ 








3.25 to 3.49 

3.50 to 3.94 

3.75 to 3.99 ___ _ 
















1 






1 


4.00 to 4.26 

4.25 to 4.49 








1 






1 


Total 








50 


148 


12 


3 


213 


Mean 

Median 

50% interval* 


1.33 

1. 27 
1. 07 to 1. 59 


1. 60 

1. 49 

1. 27 to 1. 82 


1.31 

1.30 

1. 08 to 1. 45 


1. 04 

1. 06 

0. 94 to 1. 16 


1.51 
1.43 

1. 17 to 1. 74 



*First quartile to third quartile. 



Appendix C. Evaluation and Analysis of the Unknown Parameters in the 
Life Equation with Respect to Companies and Bearing Types 



Equation (2a) of the main text, 



6.1. Summary 



Vf^D^y 



expresses the dependence of fatigue-life L on the design characteristics of the deep-groove 
bearing (Z, D a ), the bearing load P, and the "workmanship factor" f c . This appendix out- 
lines the statistical methods that were used (a) to determine "best" empirical values for the 
parameters f c , a 1} a 2 , and p of this life formula, (b) to derive the associated intervals of uncer- 
tainty, and (c) to answer various questions about the values of these parameters, from the 
basic endurance data furnished by the ASA Subcommittee, which are summarized in appen- 
dix A. These methods are applied separately in each case to the rating life L 10 and median 
life L 50 values derived from the endurance test data as described in appendix B. 

Section 6.2 summarizes the application of the statistical methods used to determine best 
empirical values and intervals of uncertainty for the parameters f c , a lf a 2 , and p in the case of 
deep-groove bearings. Sections 6.3 to 6.6 outline the statistical analyses employed to answer 
various questions about the values of f c , a x , a 2 , and p for the ball bearings of companies A, B, 
and C. In particular, section 6.3 gives the analysis employed to determine whether values 
of these four parameters are the same for the bearings of the three companies. This analysis 
is carried out separately for rating life L 10 and median life L 50 , and the postulated "complete 
between-companies homogeneity" is not supported in either instance. Section 6.4 gives the 
analysis appropriate to determining whether the data are consistent with the supposition 
that the value of p is the same for the three companies (regardless how the values of the other 
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parameters may differ); this analysis is applied to the L l0 and L 50 data, with an affirmative 
conclusion in both instances. It is concluded further that the data are consistent with the 
supposition that this common value for p is equal to 3. 

Section 6.5 presents the analysis employed to determine whether the values of the param- 
eters J c , a u a 2 , and p are the same for the three types of deep-groove bearings (B-l, B-2, and 
B-3) for which the data from company B were available. A negative conclusion is reached. 
These data are then reanalyzed to determine whether they are consistent with the supposi- 
tion that p=3 for each of the three types, regardless of differences in the other parameters. 
An affirmative conclusion is reached in this case. 

Section 6.6 concerns (a) the extent to which the L 10 and L 5Q values are consistent with the 
supposition that the values of a u a 2 , and p are 2/3, 1.8, and 3, respectively, as given in [15]; 
and (b), the determination of more precise values for f c in those cases in which the foregoing 
supposition is supported. 

6.2. Determination of "Best' 1 Values for the Parameters and Their Associated Intervals of 

Uncertainty 

As shown in section 2.2 of the main text, if natural logarithms are taken of both sides 
of the life equation (eq (2a)), the resulting equation expresses the logarithm of rating life L l0 
(or median life L 50 ) as a linear function of the logarithms of the characteristics of the bearing 
(Z,D a ), and the bearing load P, with coefficients that are simple functions of the "workmanship 
factor" / c and the exponents, a u a 2 , and p; i. e., 

Y=b +b l x l + b 2 x 2 +b 3 x 3 , (Cl) 

where 

F=ln L ~ 



£i=ln Z 
x 2 =\n D a 
x 3 =\n P j 

and 

b =p\nf c =pa 

bi=pai 
b 2 =pa 2 



(C2) 



(C3) 



are unknown constants to be estimated from the data. 

The variables x h x 2 , and x 3 are fixed variates. Their values are uniquely determined by 
the design of the bearing and the bearing loads that are used in the tests. The variable F, on the 
other hand, denotes the mean values of In i 10 , or In L b0 , for the population of all bearings with 
characteristics x t and x 2 , tested at load x 3 . 

In the practical situation Y is never known, but must be determined from the results 
of endurance tests. The methods used for obtaining such estimates of Y from endurance-test 
data are given in appendix B. To distinguish Ffrom an empirical estimate of it, the estimate 
will be denoted by the lower-case letter y. 

Generally speaking, an estimate y is a random variable, having a probability distribution 
that depends on x u x 2 , and x z . We assume that the mean of this distribution is Y= Y (x h x 2 , x 3 ), 
and that its dispersion, or more precisely, its variance, is inversely proportional to the r umber 
of bearings w in the test group from which the estimate y =y (x Xl x 2 , x; w) is derived (cf. 

305 



appendix B) ; that is, 



mean of y—Y 
variance of y=- 



4 

w J 



(C4) 



where a 2 denotes some positive cor stant, and w is the number of bearings in the test group. 

The statistical methods used to estimate the unknown parameters b , b u b 2 , and b s from the 
data are termed regression techniques. The books by Anderson and Bancroft [1], Dixon and 
Masse}' [3], Hald [10], Kempthorne [11], Mood [16], and Wilks [20] give extensive discussion 
of these techniques. For completeness, some of the techniques and rationale of regression 
analysis bearing on the work embodied in this report are summarized below. More detailed 
discussions can be found in the above references. 

Estimation. The problem of estimating the unknown parameters in the life equation can 
be stated as follows: Given independent observations (y a ; x la , x 2a , x Za ; w a ) from n test groups 
(a=l,2, . . . , n) where y a is the estimate of In L 10 or In L 50 , x ia , x 2a are the logarithms of the 
bearing characteristics Z and D a , x 3a is the logarithm of the load, and w a is the number of 
individual bearings tested for the ath test group ; required to estimate the values of the param- 
eters boj &i, b 2) and 6 3 , in eq (Cl) using some optimum method of estimation. 

Estimates for the b t (i=0,l,2,3) that are free of systematic error and have smaller variances 
than any other linear unbiased estimates are obtained by minimizing the quadratic form, 



Q=^ l w a (y a ~bo~b l x la —b 2 x 2a ~b z x Sa ) 2 

a = l 



(C5) 



with respect to each of the b t (i= 0,1,2,3). 

The resulting normal equations defining the parameter estimates can be written in the 
form, 15 



where 



aoobo+aoib 1 +ao2b 2 +aozh=go 

A A A A 

a 10 bo+anb 1 +ai2b 2 +a n h=gi 

A A A A 

a 2 obo J ra 2 ib i +a 22 b2+a 2Z b z =g 2 

A A A A 

a z obo J ra 3 ibi+a Z2 b 2 +a 3Z b z =g 3) ^ 

n 

aij=^-iW a x ia x ja i,i,=0,l,2,3- 



(C6) 



(C7) 



gi=^jW a x ia y a ^=0,1,2,3. 



*=i 



and £ a=l for all a. 

If one defines the residual sum of squares by 



n r 3 A -12 

a = l L i=0 J 



which also can be written in the alternative form 



then the optimum estimate of a 2 is 



a=l i=0 



s 

n—4c 



(C8) 



(C9) 



(CIO) 



15 A caret (a) is used here to distinguish the fact tnat the solutions of the normal equations are estimated values of the parameters and not 
the parameters themselves. 
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The quantity (n— 4) is the rank of the quadratic form S, and is termed the degrees of freedom 
associated with S. 

The above method of estimation does not depend on y a having a particular assumed 
probability distribution. All that is necessary to specify about the probability distribution of 
y a is that it possess a finite mean and variance. The solution of the normal equations (C6) 
has the property that the weighted sum of squares of the deviations about y a will be a minimum. 
This is a so-called "least-square" property of the solutions, but it is only a consequence of the 
method and is not the justification for using this method of estimation. The justification for the 
method is that this is the only one that results in minimum variance unbiased estimates for the 4 
6, (2 = 0,1,2,3). 

Values for ft, (i=0, 1, 2, 3), and hence for a t and p, can be obtained by using either the 
In Lio or In L 50 values for y a . In all cases where the a, and p have been obtained for rating 
life Lio, another set of parameters have also been calculated for median life L 5 o. 

Intervals of uncertainty and inferences. The methods so far discussed for finding estimates 
of unknown parameters need no assumption as to the form of the underlying probability dis- 
tribution of y a . However, something more must be assumed about the distribution of y a if 
(a) one wishes to place an interval about an estimate of a parameter that will include the 
"true" (or population) value of the parameter with given assurance, or (b) if one desires to 
make inferences about the parameters of the life equation for the population from which the 
bearings are a sample. 

Although the endurance lives for individual bearings may follow a Weibull distribution, 
the distribution of y a will not be of this form. However, the estimate y a (cf. appendix B) is an 
average of several independent estimates, each based on linear functions of six or less order 
statistics. Hence, by the central limit theorem, the distribution of the estimate y a will be 
approximated by a normal distribution when n is large (cf. Cramer [2, p. 213]). The statistical 
tests of significance used in this report are not greatly affected by moderate departures from 
normality. Therefore, for making all inferences, it will be further assumed that the estimates 
y a follow a normal distribution. 

The intervals of uncertainty calculated for each parameter are 95-percent confidence 
limits (equivalent to the usual "two-sigma" limits), which were referred to in the main text. 
Confidence limits for the parameter p=—b 3 can be calculated by using conventional methods. 
However the confidence limits for a , a u and a 2 are somewhat more complicated. The method 
used here, for this purpose, is sometimes referred to as Fieller's theorem (cf. Fieller [4]). 

In order to make inferences about the parameters with respect to the different companies 
or bearing types, certain statistical tests of significance were used in this paper. These are all 
based on a test statistic F, termed the variance ratio or F-ratio, which takes the form 

Jpr = QlM 

where Qi and Q 2 are quadratic forms calculated from the data and v 1 , v 2 are the respective ranks 
of the quadratic forms. The explicit expressions for Qi and Q 2 depend upon the particular 
hypothesis being tested. The subsequent sections that employ a variance-ratio statistic also 
give the explicit expressions for the two quadratic forms. 

If the hypothesis being tested is true, then the calculated variance ratio will deviate from 
unity in accordance with its tabulated distribution. However, if the hypothesis is false, then 
the variance ratio will be substantially greater than unity, and the "more false" the stated 
hypothesis, the larger the value for the variance ratio. Thus when the hypothesis tested is 
false, this will be detected by an abnormally large .F-ratio. In order to objectively determine 
whether a calculated variance ratio is significantly greater than unity, one selects from tables of 
the variance-ratio distribution a critical value of F, such that there is only a small probability 
of the calculated variance ratio exceeding the critical value from purely chance causes. The 
critical value for F used for all variance-ratio tests, in this paper, has been selected so that there 
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is only a probability of .05 of its being exceeded by a calculated variance ratio from purely 
chance causes. This critical value will be denoted by F^{vi,v 2 ). More extensive discussions of 
these procedures are given in Kempthorne [11, chap. 5]. 

All statistical tests of significance to test relevant hypotheses have been carried out both 
for rating life Zio and median life L 50 . 

6.3. Analysis to Determine Whether Companies Have Common Values for All the Parameters 

in the Life Equation 

This section deals with the details of estimating the values of the parameters in the life 
equation for each company. Furthermore, a statistical analysis is made to determine whether 
the companies have common values for all the parameters in the life equation. The 7^-ratio 
(eq (Cl6)), which is used to test this hypothesis, is obtained from the following procedures: 
A single set of parameters, b t (i=0, 1,2,3), is obtained by first fitting all the data, irrespective 
of company, to the logarithmic life equation, and then calculating the resulting residual sums 
of squares S (eq (C14)) having 206 degrees of freedom. If the hypothesis of common values 
for all the parameters is not true, then a better fit to the data can be made by fitting the life 
equation separately to each company. These calculations result in the individual residual sums 
of squares Si, S 2 , and S 3 (eq (C13)) having 46, 144, and 8 degrees of freedom, respectively. 
Thus the total residual sums of squares (S1+S2 + S3) will have 46 + 144+8 = 198 degrees of 
freedom. Then the difference between S and (Si+#2+&), i. e., {S— ($i+£ 2 +S 3 )}, is also a 
quadratic form having 206 — 198 = 8 degrees of freedom. If a substantially better fit was ob- 
tained by fitting a separate life equation to the data for each company, as compared to a single 
life equation, the difference between the two residual sums of squares {S— (S1+/S2+S3) } will be 
large. To determine whether this difference is statistically significant the variance ratio 
(C16) is employed. 

Mathematical formulation. It will be convenient to adopt the following notation: Let the 
superscript (u) = l, 2, 3 refer to companies A, B, and C, respectively. Also, for each company, 

A 

let b$ u) C£=0,l,2,3) refer to the estimates of the parameters in eq (Cl), and a\ u >, g\ u) denote 
the sums of cross products defined in eq (C7) and (C8) 16 . Then the normal equations that give 
the parameter estimates for the wth company are 

J2a^h u) =gl u) i=0,l,2,3, (Oil) 

i=o 

and the estimates for the parameters a t (i=0, 1, 2) and p, in the life equation, are obtained from 
the relations given by eq (C3). These results are summarized in tables 6 and 3, respectively, in 
the main text. 

The hypothesis that all parameters in the life equation are the same for each company is 
equivalent to the hypothesis that 

ft<i) == ft(2) == jj8) i=o,l,2,3. (C12) 

Define the residual sum of squares for the uih company by 



and let 



S u =T;wuayla-?b M u) gl u) u=l,2,3, (C13) 

a=l 1 = 

S=SE«L-Sk ffi=S^" > , (Ci4) 

U = l a=l 1=0 U = \ 



where b t (i=0,l,2,3) are the estimates obtained from the solutions of the normal equations 



16 As it is only the ratio of the numbers of observations that is important for weighting, the weights w used in the calculation of the sums of 
cross products in all the analyses have been taken as integral multiples of 5; e. g., if the number of bearings in a test group is 26 (say), then w=5. 
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without regard to company differences, i. e., 

S |~E <>1 &,=£ ,/{»> t=0,l,2,3. 

j=0 Lw = l J w=l 

Then the variance ratio for testing the hypothesis given by eq (C15) is 

7^ (S — Si — So — Sz)/8 

(Si+S 2 +S 8 )/198 ' 



(C15) 



(C16) 



and the critical F value is F, 05 (8, 198) = 1.98. 

Table C-l summarizes the values for S u (w= 1,2,3) and the other quantities needed to 
evaluate eq (C16). 

Table C— 1. Synopsis of calculations to determine whether companies have common values for all 

parameters in life equation 





Degrees 

of 
freedom 


L\o 


Lso 




Mean 




Mean 


Companies combined 




S u 


square 


S u 


square 


206 


446. 684633 




348. 285324 




A 

B 

C 


46 
144 

8 


102. 552212 

285. 121962 

6. 160858 




89. 734079 

201. 935797 

8. 357822 




Sum 


198 


393. 835032 


1. 9891 


300. 027698 


1. 5153 


Difference 


8 


52. 849601 


6. 6062 


48. 257626 


6. 0322 


Fromeq (C16) 




6. 6062 
1. 9891 


-3. 32 


6.0322 
F -1.5153~ 3 - 98 



The values of the /'"-ratio (eq (C16)) calculated from the quantities in table C-l are 

L 10 : F=3.32 

L 50 : 7^=3.98 

These calculated F values are both larger than the critical value, F. 05 (8, 198) = 1.98; actually 
the probability is less than .0005 of having an i^-ratio as large as those above from purely 
chance causes. Thus from the above statistical tests of significance, the conclusion is reached 
that the three companies do not have common values for all of the parameters in the life equation. 

6.4. Analysis to Determine Whether Companies Have a Common Value for the Exponent p 

The previous analysis resulted in the conclusion that the parameters in the life equation 
are different for each company. However, this does not exclude the possibility that all com- 
panies may have a common p, even though the a { (i=0,l,2) differ from company to company. 
This section discusses the analysis made to determine whether the companies have a common 
value for the exponent p. The analysis given here consisted of the following procedure. First, 
the logarithmic life equation (C17), having a common value of p, but allowing the a t to vary for 
each company, was fitted to all the data, and the resulting residual sum of squares S' (eq 
(C20)), having 200 degrees of freedom, was calculated. The total residual sum of squares 
from fitting the life equation separately to each company (allowing p to vary in addition to the 
a t ) is given by {S l +S 2 J rS^) having 198 degrees of freedom (cf. section 6.3 of this appendix). 
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Then the reduction in the residual sum of squares achieved by using a different exponent p for 
each company is{S'— (Si+S 2 +S d ) } having 200 — 198 = 2 degrees of freedom. To test whether 
this reduction in the residual sums of squares is statistically significant, the variance ratio 
(C21) is employed. 

Mathematical formulation. The logarithmic life equation, having a common value for 
the exponent p, can be written for the ath test group in the uth company as 

Y ua =b ( u) + b[ u) x[^ + ¥ 2 u) x^ + b 3 x^ a=l, 2,...,n u ; u=l, 2, 3. (Cl7) 

Note that although each company has the same parameter b 3 in the above equation, the para- 
meters b ( u) , b[ u) , and b ( 2 u) are different for each of the three companies. Thus there [are 10 
different parameters, i. e., bl u) (i=0,l,2; u= 1,2,3) and b z , to be estimated from the data. 
The normal equations for estimating these parameters are 



S aflbf+aft h=gl u) i=0, 1, 2;u=l, 2, 3 



lb1baWbr+a 3 X=fh, 
where 



(ci8) 



u=l u=l 

Thus the set of equations given by (C18) is a system of 10 linear equations in 10 unknowns. 
Once the solutions are obtained, the estimates for a ? -(/i=0,l,2) in the life equation are calculated 
from the relationship 

b { - u) 
al u) =— -ir- i=0,l,2;w=l,2,3. (C19) 

h 

The residual sums of squares (denoted by S') associated with fitting the life formula (Cl7) to 
the data is given by 

3 Tin 32 A 

S'=T, S w Ma j/L-S S 6™> gP-hg,. (C20) 

il = l a = l U=l 1 = 

Then, to test the hypothesis that the companies have a common value of p, regardless of 
the values for the other parameters in the life equation, the variance ratio 

-p_ (S' — Si— S 2 — S3) /2 /non 

*~ (S l + S t +S t )/198 ' { ^ Z1) 

having 2 and 198 degrees of freedom is used. The critical F value is F. 05 (2, 198) = 3.04. 
The values of S' for both L l0 and L 50 are 

L 10 : £"=393.272847 d. f. = 200 

(C22) 
L m : S'=301. 687871 d. f.=200, 

and the calculated 17 F values (using eq (C21)) give 

Z 10 : F=-0.14l] 

Y (C23) 

L- M : F=0.548. J 



17 From theory, the calculated value for the .P-ratios can never be negative. The reason for the negative value of F for Zio is that the value 
for the numerator of eq (C21) is only accurate numerically to one decimal place on account of round-off errors arising from the solution of the 
normal equations (C18). Thus, if the hypothesis of a common p value is true, then the F-ratio will not be large and round-off errors may affect 
the resulting calculation. Alternatively, if the null hypothesis is false, then the calculated F-ratio will be larger then 3.04 and the round-off error 
should be of no consequence. 
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Because both variance ratios are smaller than the critical value F. 05 (2, 198) = 3.04, the 
conclusion drawn from this statistical analysis is that the data support the hypothesis of a 
common value of p for the three companies. This holds both for rating life (Z l0 ) and median 
life (L 5 o). The values for the common p are given in table 2 of the main text. The values 
for the remaining constants in the life equation di(i=0,l,2) are found from the relations (( J19). 
These results are summarized in table 5 of the main text. 

6.5. Analysis to Determine Whether the Three Bearing Types From Company B Have Con- 
sistent Parameter Values 

The analyses, discussed in previous sections, dealt with determining whether there are 
differences in the parameters of the life equation between companies. This section investigates 
(a) whether three different bearing types made by company B have common values for all the 
parameters in the life equation, and (b) whether the exponents p calculated for each bearing 
type are consistent with the value of p = 3. 

The analysis for (a) is similar to the analysis made in section 6.3 of this appendix; i. e., 
separate life equations were fitted to each bearing type and the resulting residual sum of 
squares was compared with the residual sum of squares arising from fitting a single equation 
to all data from company B, irrespective of bearing type. The variance ratio for statistically 
testing (a) is given by eq (C27). 

The analysis for (b) was governed by the following considerations. If the true (or popula- 
tion) value of the exponent p is p= 3, regardless of bearing type, then the estimates for p 
obtained by fitting a separate life equation to each bearing type should not differ from p=3 by 
more than the dispersion inherent in the endurance lives of the bearings. The agreement of 
the values of p estimated for each bearing type with p=3 is tested for statistical significance 
by the variance ratio (C28). 

Mathematical formulation. The 148 test groups from company B can be divided into 
three bearing types corresponding to 37 groups for B-l, 94 groups for B-2, and 17 groups for 
B-3-type bearings. Let these types be denoted by 0=1, 2, 3, respectively. Also define 



**-ij — ^ 1 W va %ia %ja ^j J — ^, 1? A 3,' V — 1, 2, O 



fll^S^iilV,'. 1 i=0', 1, 2, 3; »=1, 2, 3, 

a=l 



(C24) 



where n v is the number of test groups for bearing tvpe v. Then the normal equations, which 
determine the estimates for the parameters in the logarithmic life equation, are 



2}A<W> = <?{•> i=0, 1, 2, 3; v=l, 2, 3. 



(C25) 



Thus for each bearing type, the values for the parameters a*(i = 0,l,2) and p in the life equation 
can be estimated from the relations 



b iv) 

PfJV) — - - l 

a i a 



£(*>=_ J<*>. 
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Therefore the residual sum of squares for the pth bearing type is 



nv 3 A 



(C26) 



having (n„— 4) degrees of freedom. Then the hypothesis of common parameters for the three 
bearing types can be tested by the variance ratio, 



F= 



(£ 2 -£ (1) -S (2) -£ (3) )/8 
(S (1) 4-S (2) + S (3) )/136 ' 



(C27) 



having a critical value of ^.05(8, 136) =2.01. 

Since the analysis given in section 6.4 of this appendix reached the conclusion that all 
companies have a common value for p, and since this value (cf. table 2 of main text) with its 
associated uncertainty includes the value p = 3 given in [15], it seems also desirable to test a 
second hypothesis that the value of p for each bearing type is consistent with p=3. The 
F-ratio for this hypothesis is given by 



ID 



fp (1) -3) 2 , (ft (2) -3) 2 , (|> (3) -3)' 



no 

^3J 



] 



(S (1) + S (2) + S (3) )/136 



(C28) 



33 



where the p iv) (v=l,2,3) refer to the estimates of p obtained for each bearing type, and 
(#=1,2,3) is the element occurring in the last row, last column of the inverse matrix to WAffW 
(t>=l,2,3). The critical value for the variance ratio (C28) is F M (3, 136) =2.67. 

The estimates of a\ v) (i=0, 1, 2) and p (v) obtained from the solutions of the normal 
equations (C25) are summarized in tables 7 and 4, respectively, in the main text. The calcu- 
lations for the variance ratio (C27) are summarized in table C-2. 

Corresponding to the hypothesis that the three bearing types have the same parameters 
in the life equation, the calculated variance ratios (C27) yield 



(C29) 




Table C-2. — Synopsis of calculations to determine whether all bearing types have common parameters in life 

equation (company B only) 






£10 


Lm 


Type of bearing 


Degrees 

of 
freedom 






£<•> 


Mean 


£<•> 


Mean 








square 




square 


All types combined 


*144 


*285. 121962 




*201. 935797 




B-1 


33 


68. 558298 




33. 016973 




B-2 


90 


174. 325939 




131. 142017 




B-3 


13 


11. 537480 




8. 505638 




Sum 


136 


254. 421716 


1. 8707 


172. 664628 


1. 2696 


Difference 


8 


30. 700246 


3. 8375 


29. 271169 


3. 6589 


Fromeq (C27) 


'-!SE-™ 


3.6589 
F= T2696 =2 - 88 



*From table C-l. 
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Because the critical lvalue is F^ b (8, 136) = 2.01, both the L ]0 arid L 50 calculated variance 
ratios are statistically significant. Therefore, one could conclude from the above F-ratios that 
the data support the hypothesis that parameters do differ between bearing types. However, this 
does not exclude the possibility that the values for the exponent p are consistent with the value 
p = S. Substituting the appropriate quantities in eq (C28) results in 

L 10 : F=1.18*\ 

r (C30) 

L 50 : F=2.32J , 

where 

Off =0.53294 

QV =0.107395 

Qf =0.396744. 

Thus, as both of the above calculated F-ratios are less than the critical value, F.05 (3,136) 
= 2.67, the conclusion can be made that the values oj p are consistent with p=3 for different 
bearing types, although possibly some (or all) of the other parameters in the life formula (i. e., 
a*(i=0,l,2) may differ among the different bearing types. 

The values for a\ v) (i=0,l,2; v= 1,2,3) arising from the analysis by ball-bearing types 
have very large confidence limits (intervals of uncertainty). This is mainly due to the fact that 
an analysis restricted to one bearing type is essentially an analysis on bearings having almost 
the same values for Z and D a . In order to estimate the a\ v) with good precision, it is necessary 
to have results for bearings having wide variations with respect to Z and D a . Thus the esti- 
mates for a t based on all bearing types for company B (table 5 or table 6) have substantially 
smaller confidence intervals as compared to the intervals based only on a single bearing type. 

6.6. Determination and Analysis of f c Based on the Parameter Values 
a 1 =2/3 / a 2 =1.8, andp=3 

The values for the parameters a 1 , a 2 , and p given in [15] are ^=2/3, a 2 =1.8, andp=3. 
If these parameter values are valid for the data at hand, then more precise estimates for the 
" workmanship " parameter a (or/ c ) can be made for each company or bearing type. These 
will generally have better precision compared to the estimates of a made when the other 
parameters in the life equation are simultaneously estimated along with a . This section 
considers the problem of verifying whether the parameter values, given above, are valid for the 
given data, and for those cases where this is true, estimates of a (or/ c ) are obtained assuming 
these values for the other parameters. 

The procedure for determining whether the values a! = 2/3,a 2 =1.8, and p=3 are valid for a 
given classification of the data (with respect to a company or bearing type) is to fit the data to 
the life equation, using the assumed values for a h a 2 , and p. Thus there is only one unknown 
parameter, a , in the life equation to be estimated. Then the resultant residual sum of squares, 
denoted by R (eq (C37)), can be calculated having (n—1) degrees of freedom. Alternatively, 
the life equation can be fitted to the data such that all the unknown parameters are simultane- 
ously estimated. The residual sum of squares from this latter fit, S, will have (n — 4) degrees 
of freedom. Then, if the above parameter values are not consistent with the given data, R 
will be appreciably larger than S. The variance ratio (C38) is used to determine if the differ- 
ence between these two residual sums of squares, i. e., (R—S) having [(n—l) — (n—4)] = 3 
degrees of freedom, is statistically significant. 

Mathematical formulation. Let n be the number of test groups within a particular classifi- 
cation (either by company or bearing type). Then, assuming the values a 1= =2/3, a 2 =l.8, and 
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p=3, the logarithmic life equation for the aih. test group can be written as 

Y a =3{a J rx a } a=l,2, . . ,,n } 

where 



The resulting normal equation for estimating a Q is 



(C31) 



(C32) 



do= 



_a = l O « = 1 



(C33) 



which also can be written as a function of the sums of cross products g i} a ih e. g., 



aoo \ 3 v 



aoi + 1-8 a 02 — a 03 



.)}■ 



Hence the estimated variance of a is 



variance (oo)=^ — f 
9a o 



where 



R 



9 S ^ (y — #«) -a^al 



(C34) 



(C35) 



(C36) 



n — 1 



n-1 



The residual sum of squares R, having (n—1) degrees of freedom, can also be written as a func- 
tion of the sums of cross products, 



i? = 9 



{1 n 4 2 4 4 a "1 

- ]C w a yi—- #i — 1.2 # 2 +- # 3 +- an+3.24 022+^33+ 2.4 ai 2 — ^ ai 3 — 3.6 a 23 — floo^o I . 



(C37) 



In the analyses made in the preceding sections, the ball-bearing data have been analyzed 
with respect to individual companies or bearing types. It thus seems desirable to determine 
whether the data within these classifications support the hypothesis that a x = 2/3, & 2 =1.8, and 
j9— 3. The variance ratio used to test this hypothesis is 



S/(n-4)' 
where, for testing within companies, 

' Si for company A 

S=< S 2 for company B 

h Si for company C _, 

and for testing within bearing types, af] is replaced by Af] , 

C 8 a) for B-l 



(CSS) 



>■ defined by eq (Cl3) and given in table C-l, 



S=-l 



S i2) for B-2 
L S (3) for B-3 J 



defined by eq (C26) and given in table C-2, 



314 



and R (eq (C37)) refers to the calculated residual sum of squares within the particular 
classification. 

Table C-3 summarizes with respect to companies A, B, and C the values of a , 11, the 
calculated F-ratio (eq C38), and the critical F value. Table C-4 summarizes the same quan- 
tities for the B-l, B-2, and B-3-type bearings made by company B. 

The results of this analysis show that the values ai=2/3, a 2 =l.S, and p=3 are consistent 
for rating life L 10 data with respect to each of the three companies. However, for median 
life L 50 , these assumed parameter values are consistent only for the data from company A. 

The results for the same analysis made on the three different bearing types, indicate 
that the assumed parameter values are consistent only for the B-l and B-2-type bearings 
witli respect to median life Z 50 . 

It is interesting to note that the analysis for company B (ignoring bearing types) showed 
that the assumed values for the parameters are consistent with the L 10 data. However, a 
finer analysis by bearing type revealed that these values are not valid for the B-3-type bear- 
ings. This apparent inconsistency stems from the fact that the analysis for company B, 
taken as a whole, is dominated by those bearing types having the larger number of test 
groups, i. e., B-l and B-2, and for these types the parameter values were found to be 
consistent with the Z 10 data. 

The estimates for f C) assuming ai=2/3, & 2 =1.8, and p = 3, are summarized in table 8 of 
the main text for rating life L i0 . This summary also includes the value of j c for company 
D computed from only three test groups. (Because of the small number of test groups, it 
was not possible to verify whether the assumed parameter values are valid for these data.) 



Table C-3. Summary of computation* (by companies) to test hypothesis that data are consistent with assumed 

values 0,1 = 2/3 } a 2 = 1.8, p = 3 



Company 


Lio 


L 50 


Critical 
F 


&o = h\f c 


R 


F 


d = \nfc 


R 


F 


A__ 


8. 4205 
8. 5021 
8. 1001 


102. 539286 

287. 613198 

8. 375562 


*-0. 06 

0. 42 
0.96 


8. 9382 
8. 925 1 
8. 5832 


94. 367916 
217. 560672 

8, 488584 


0.79 
3. 71 

9. 92 


2. 80 
2. 65 
3.59 


B 

C 



* Negative value due to rounding in calculations. 



Table C-4. Summary of computations (by bearing type) for company B to test hypothesis that data are consistent 

with assumed values a\ = 2/3, a 2 = 1.8, p=3 



Type 


Lio 


Lw 


Critical 
F 


do=ln/ c 


R 


F 


d = lnf c 


R 


F 


B-l 

B-2 


8. 4575 
8. 5236 
8. 5203 


73. 821978 

181. 128312 

27. 014247 


0. 84 

1. 17 

5.81 


8. 8636 

8. 9482 

9. 0022 


44. 883126 

141. 731604 

18. 905517 


3.95 
2. 42 
5. 21 


2. 8G 

2. 70 

3. 24 


B-3 
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